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G\ ', ABSTRACT 

G\ ■ 

i \ We examine the relation between presupernova stellar structure and the distribution 

| of ejecta in core-collapse supernovae of types lb, Ic and II, under the approximations 

■^■j- ■ of adiabatic, spherically symmetric flow. We develop a simple yet accurate analytical 

formula for the velocity of the initial forward shock that traverses the stellar envelope. 
j> ■ For material that does not later experience a strong reverse shock, the entropy deposited 

by this forward shock persists into the final, freely-expanding state. We demonstrate 
that the final density distribution can be approximated with simple models for the 
final pressure distribution, in a way that matches the results of simulations. Our results 
indicate that the distribution of density and radiation pressure in a star's ejecta depends 
on whether the outer envelope is radiative or convective, and if convective, on the 
composition structure of the star. 

Our models are most accurate for the high-velocity ejecta cast away from the pe- 
riphery of a star. For stellar structures that limit to a common form in this region, the 
resulting ejecta limit to a common distribution at high velocities because the blast wave 
forgets its history as it approaches the stellar surface. We present formulae for the final 
. | density distribution of this material as a function of mass, for both radiative and effi- 

ciently convective envelopes. These formulae limit to the well-known planar, self-similar 
solutions for mass shells approaching the stellar surface. However, the assumption of 
adiabatic flow breaks down for shells of low optical depth, so this planar limit need not 
be attained. The event of shock emergence, which limits adiabatic flow, also produces 
a soft X-ray burst of radiation. Formulae are given for the observable properties of 
this burst and their dependence on the parameters of the explosion. Motivated by the 
relativistic expansion recently inferred by Kulkarni et al. (1998) for the synchrotron 
shell around SN1998bw, we estimate the criterion for relativistic mass ejection and the 
rest mass of relativistic ejecta. 

We base our models for the entire ejecta distribution on the high-velocity solution, 
on our shock velocity formula, and on realistic radiation pressure distributions. We also 
present simpler, but less flexible, analytical approximations for ejecta distributions. We 
survey the ejecta of the polytropic hydrogen envelopes of red supergiants. Our models 
will be useful for studies of the light curves and circumstellar or interstellar interactions 
of core-collapse supernovae, and the birth of pulsar nebulae in their ejecta. 
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Subject headings: gamma rays: bursts — hydrodynamics — shock waves — supernovae: 
general — supernova remnants 



1. Introduction 

All stars of initial mass greater than about 8M end violently when their cores collapse; 
many or all of these produce supernova explosions. Supernovae are of paramount importance to 
the evolution of galaxies: they form new elements and recycle the products of stellar evolution. 
Their energy is converted into motion and heat, affecting the dynamical state of the surrounding 
gas; indeed, they are responsible for much of the hot, X-ray emitting gas in galaxies, and for the 
formation of superbubbles and outflow in starburst galaxies. Observationally, supernovae provide a 
natural laboratory in which to test theories of stellar evolution and of the circumstellar environments 
of massive stars: supernova 1987A is a notable example ( McCray 1993| ). 



How does the stellar structure of a progenitor affect the aftermath of its explosion? The 
observable properties of a core-collapse supernova are determined by the distribution of material 
and temperature in its freely-expanding ejecta. The supernova's optical display begins with a 
photospheric stage, in which a recombination front moves inward through the expanding ejecta. 
This process is pre-determined by the temperature, density and opacity profiles of the outcast 
material. At some point, the energy from radioactive decays begins to dominate the light curve. 
The nebular display is determined by the distributions of both decaying isotopes and gamma- 
ray opacity. Later, as the ejecta collide with circumstellar or interstellar material, two shocks 
are born at the interface: a forward shock accelerating the surrounding material, and a reverse 
shock decelerating the ejecta ( [VlcKee 1974| , Chevalier 1982 , Dwarkadas fc Chevalier 1997 , Truelove 



& McKee 1998). This interaction depends critically on the distributions of ejected and ambient 
material. It is responsible for the appearance of the supernova in the radio, X-ray and infrared 
bands for centuries after the explosion. 

The distributions that characterize a supernova's ejecta are, in turn, determined uniquely by 
the structure of the progenitor star (although the distribution of radioactive elements is an added 
complication). The relation between progenitor and ejecta is established by the hydrodynamics of 
the supernova explosion itself, which occurs in three stages. In the first stage, a shock propagates 
through the stellar material to be ejected, pressurizing it and setting it into motion: we call this 
the blastwave stage, following Chevalier (1976). Once the shock arrives at the stellar surface, a 
rarefaction phase begins: the material accelerates as its heat is converted into kinetic energy. This 
stage (the most difficult to model) ends when acceleration ceases, once the ejecta's pressure is 
too slight to be dynamically important. In the ejecta stage, the outcast material was coasts into 
homologous expansion, with a velocity distribution v = r/t. 



The correspondence between progenitor and ejecta can be established numerically on a case- 
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by-case basis. This approach has been pursued for many years (e.g., Falk & Arnett 1977). However, 
many aspects of the late stages of stellar evolution are uncertain, including rotation and pulsations 
(Heger et al. 1997), nuclear reaction rates, mass loss prescriptions and models for semiconvection 



( Woosley, Langer fc Weaver 199~3| ), and it is likely that interactions with stellar (e.g., Nomoto, 



Iwamoto, & Suzuki 1995) or substellar companions will affect supernova progenitors. Numerical 
simulations may describe specific cases quite accurately, but only limited insight into the physical 
relationship between progenitor and ejecta can be derived from an individual case. 

Supernova progenitors and their ejecta display strikingly similar features, when the initial and 
final densities are considered functions of a Lagrangian variable, like the enclosed mass. Does this 
imply that there is a simple relationship between the two? This paper will investigate the existence 
of such a relationship, that is, the degree to which the variables that describe a mass shell in the 
progenitor control its properties after the explosion happens. A local relation between the initial 
and final variables implies that the physics of the explosion can be understood in simple terms, 
despite the complicated hydrodynamical interactions involved. The pressure-based models of §|6| 
demonstrate that this is the case. 

Another goal of this paper will be to provide simple models for broad classes of supernova 
progenitors, like red supergiants. When progenitor stars are similar in structure, their ejecta will 
also take similar distributions. A simple model that can reproduce these ejecta accurately enough 
to show the differences between them, like the harmonic-mean models of §|7[ should prove valuable 
for studies of supernova light curves and circumstellar interactions. 



1.1. Previous investigations 

The analytical study of supernova hydrodynamics has relied on familiar self-similar analogs. 
Over a considerable fraction of the progenitor, the density distribution resembles a power law in 
radius: po oc r^ n with n ~ 1.5 — 3 above the core, but usually with a shallower distribution 
further out. In a strictly power law, spherical distribution, the blastwave stage takes the self- 
similar form given by Taylor (1950) and Sedov (1959). Chevalier (1976) likened the structure of 
a blastwave in an increasingly steep density distribution to a progression of Sedov blastwaves in 
different power law media, and concluded that a detached shell (density maximum) would occur if 
density became steeper than a critical value. By analogy to the Sedov blastwave with po oc r^ 3 , 
Bethe & Pizzochero (1990) conjectured that the shock in supernova 1987A would have nearly a 
constant velocity. Chevalier &: Soker (1989) approximated the central progenitor distribution as 

17/7 

po oc r , which produces a particularly simple blastwave distribution (the Primakoff-Sedov 
form) in which density and pressure are power laws in radius and time. 

The resemblance to a Sedov blastwave must break down as the shock approaches and emerges 
from the stellar surface, and the rarefaction phase begins. But there is an analytical solution to 
guide the understanding of this process as well. The density varies as a power of depth in the outer 
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layers of a star, and radial curvature can be neglected near enough to the stellar surface. A planar 
version of this problem was initially considered by Gandel'man & Frank-Kamenetsky (1956) and 
solved by Sakurai (1960). In this context, the shock front accelerates as it runs down the declining 
density gradient, and material expands into vacuum after the shock emerges from the surface. The 
ratio of final to postshock velocity is the same for all fluid elements in this solution; the density and 
pressure approach steep power laws in velocity as the material expands freely. This solution has 
been applied to the pulse of emission associated with shock breakout and the early phase of the light 
curve dChevalier 1992 ), as well as the distribution of high-velocity ejecta from the explosion (e.g., 



Chevalier fc Soker 198S , Imshennik fc Nadezhin 1989] , Litvinova fc Nadezhin 199C| ). Imshennik & 



Nadezhin (1989) and Chevalier (1992) showed that spherical curvature should complicate matters, 
because a mass shell in the planar solution must travel a long way (compared to its initial depth) 
to reach its final velocity. Kazhdan & MurzinaQ (1992) took the dynamical effect of sphericity into 
account in a first-order perturbation analysis of Sakurai's solution. In this paper, we will find the 
correction of the planar solution to higher order by a different method. 

Sakurai's planar solution describes the end of shock motion and the beginning of rarefaction. 
But even with the first-order effect of sphericity accounted for, this solution only holds for the 
highest-velocity ejecta, or for the brief period in which the shock crosses and ejects a thin outer 
shell. It does not describe how this brief period joins onto the blastwave phase that preceded it and 
the rarefaction phase that follows it, in which the bulk of the ejecta approaches its final velocity. 
The rarefaction phase is very complicated in general, even if the blastwave was one of the Sedov 
forms; rarefaction is therefore a formidable problem. 

There is one counter-example. For the special case of the Primakoff blastwave, in which density 
and sound speed are power laws of radius and time, a rarefaction wave will travel along a power 
law trajectory. Chevalier & Soker (1989), using the Primakoff form to approximate the blastwave 
in SN1987A, argued that the spherical rarefaction process will be self-similar in this case. Then, 
the ratio between the final velocity and the velocity at the rarefaction wave should be the same 
for any mass shell. This simple argument specifies that the density should settle into a shallow 
power law of velocity, in the inner region of the ejecta. Joining this shallow inner power law with 
the steep outer power law specified by Sakurai's planar solution, and imposing overall mass and 
energy conservation, Chevalier & Soker found the velocity of the break between the two behaviors. 
Although it is based on idealizations of the flow, this model compares admirably with simulations 
involving blue supergiant progenitors. 



1 We are grateful to the referee, Roger Chevalier, for bringing this work to our attention. 
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1.2. This work 

In this paper we address each stage of the expulsion of stellar envelopes in core collapse super- 
novae, in order to arrive at a reliable model for the final ejecta distribution that can accommodate 
possible variations in progenitors' stellar structures. Wherever possible, we make use of scaling 
relations to simplify the problem. When considering red supergiants we concentrate primarily on 
the outermost composition layer (the hydrogen envelope) where our models will be most accurate, 
approximating these envelopes as polytropes when it is useful to do so. Table |l| gives an overview 
of the models developed in this paper; these are motivated by and then checked against simulations 
for which the initial configuration is the result of a stellar evolution study. These progenitors were 
kindly supplied to us by Stan Woosley and Ken'ichi Nomoto. 

We begin by examining in §|| the scaling relations that connect the hydro dynamical behavior 
of similar initial structures, under the idealization of adiabatic flow. The scaling to similar solutions 
will be especially useful in §||, where we present solutions for the high-velocity ejecta of convective 
and radiative envelopes with a single form for each. 

Section |3] discusses the general features of supernova progenitor stellar structures. These are 
marked by composition zones of vastly different density and radial scale. The blastwave shock 
typically decelerates inside each layer, and accelerates in the steep density declines between zones. 
The shock wave accelerates as it approaches the surface of the star, at the end of the blastwave 
phase. 

The models we develop will require an approximation for the detailed motion of the blastwave 
shock front. The forward shock defines the outer boundary of the flow in the blastwave stage. 
Moreover, the blastwave shock velocity sets the entropy distribution for the ensuing adiabatic flow 
(for material not subject to reverse shocks). The shock velocity thereby affects the rarefaction phase 
as well, and the relation between pressure and density in the final ejecta phase. The shock velocity 
is also necessary to set the coefficients of the power laws that describe the highest-velocity ejecta. 
In reality, the shock velocity must be determined self-consistently by matching it to the postshock 
pressure in the blastwave phase. Previous studies (Kompaneets 1960, Laumbach Sz Probstein 
1969, Whitham 1974, Klimishin & Gnatyk 1982, Koo & McKee 1990) have suggested a number 
of formulae for the velocities of shocks in inhomogeneous, three-dimensional distributions; all of 
these formulae have drawbacks in the context of supernovae, and we will suggest another simple 
shock velocity formula better adapted to this context in §|j. This shock velocity approximation is 
the product of well-known scaling laws for spherical decelerating shocks and planar accelerating 
shocks. Despite its simplicity, this approximation tracks the successive periods of shock acceleration 
and deceleration with remarkable accuracy. This approximation gives the entropy distribution in 
the ejecta, except where a reverse shock interferes. It also gives the coefficients that describe the 
highest-velocity ejecta. 

In §|5| we address the dynamics of material near the surface in the progenitor, the material 



that is ejected at high velocity. The initial distribution of this material, described in §3.1, has a 
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predictable form in radiative and polytropic envelopes. Having a reliable estimate of the velocity 



of the shock that launches this material, we calculate in §5.1 the coefficients of the power laws that 
describe its ejecta, as given by the planar, self-similar solution (Gandel'man & Frank-Kamenetsky 
1956, Sakurai 1960; see the Appendix). Most material attains a smaller velocity than the planar 
formula would give, because spherical expansion reduces the pressure below its planar value. This 
is a geometrical effect, however, and because the initial density distributions of this material are 
similar among broad classes of envelopes (radiative or convective), the final ejecta distributions are 
also similar. In § |5.2| , we present simple formulae that capture the deviation of important quantities 
away from their self-similar power laws. 

The self-similar solutions that give power laws for the high-velocity ejecta are not valid to 
infinite velocity, of course. The assumption of adiabatic flow breaks down for material involved in 
shock breakout, because radiation diffuses or streams out of the gas. This process also produces a 



burst of soft X-rays. In §5.3.1 we present formulae for the upper velocity limit for adiabatic flow. 
We use our non-relativistic analysis to derive an approximate criterion for relativistic mass ejection 
and the rest mass of relativistic ejecta (if any) in § 5.3.2| . In § |5,3.3 , we present formulae for the 



characteristics of the X-ray burst. Although previous authors have treated the radiation outburst 
in a similar manner (e.g., Imshennik fc Nadezhin 1989 ), the shock velocity formula of §|| allows us 



to elucidate the effects of stellar structure on the outburst. 

The results of §|5| apply to the high- velocity ejecta, those at velocities a few times greater than 
the systemic velocity. However, little of the mass or energy in the ejecta lies in this velocity range. 
For the bulk of ejecta, the details of the progenitor's structure break the scaling relations that allow 
the high-velocity ejecta to be given accurately by simple formulae. We note that the final pressure 
distribution is always smooth and continuous in the ejecta, even if the density distribution is not. 
For this reason, we base our models for the entire ejecta distribution, discussed in §|6|, on models for 
the pressure distribution that mimic the smooth pressure functions seen in simulations. We find a 



simple way to extend the pressure distributions derived in §5.2 inward in mass. The final density 
distribution is specified by the combination of this pressure model with the entropy left behind 
in the shock. These pressure-based models constitute a local theory, in that they relate the final 
density of a mass shell with its initial state. These models fail in regions that are traversed by a 
strong reverse shock, because we do not account for the entropy it deposits; this failure is worst in 
the inner composition layers of red supergiants. Even so, there is a jump in density between the 
ejected helium mantle and the ejected hydrogen envelope that is not affected by the reverse shock; 
a simple analysis of this phenomenon is presented in § |6.1| . 

Section [?] presents a different type of model for ejecta distributions, in which the shape of 
the final density is described as a simple, harmonic-mean interpolation between two power laws 
of velocity. This is a global model, since it gives an ejecta distribution (in terms of three free 
parameters) without referring to the detailed structure of the progenitor. For blue supergiants, 
we adopt the theory of Chevalier & Soker (1989) to motivate the choice of parameters. For red 
supergiants, we find that the harmonic-mean model can fit the ejected hydrogen quite accurately, 
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but it does not fit the ejected mantle because of the density jump between the two regions. 

In §|8] we develop simplified polytropic models for red supergiant progenitors of type II su- 
pernovae. We argue that these can be adequately described by two parameters that are ratios of 
masses in the progenitor. Moreover, these mass ratios are tightly correlated, so that there is essen- 
tially a one-parameter family of red supergiant progenitors. We investigate this family of models 
numerically to find the distribution of the free parameter used in the pressure-based model of §|6| 
and the three free parameters in the harmonic-mean models of §0; results are given in Table § 

Section [9| contains a comparison between our models and the ejecta of simulated supernovae 
whose progenitors are the end products of stellar evolution calculations. This comparison allows a 
tabulation of the errors introduced in the pressure-based and harmonic-mean models (Table |6|) in a 
realistic context. The models developed in this paper help to explain many of the features seen in 
numerical simulations of supernovae. This is true even for features are not reproduced accurately 
by the models, those caused by reverse shocks or by variations in the polytropic index at the outside 
of the progenitor. These features are discussed in §^ and in the Conclusions. 

In the Appendix, we develop a fully Lagrangian treatment of the self-similar shock accelera- 
tion problem solved by Sakurai. This problem has previously been addressed in two steps, with an 
Eulerian solution prior to emergence matched to a Lagrangian description of the subsequent expan- 
sion. The fully Lagrangian description facilitates analyses of the highest-velocity material. We use 



it to compare our results of §5^ to the suggestion of Litvinova & Nadezhin (1990) that spherical 
expansion essentially truncates planar, self-similar acceleration, and we find that this happens at 
about three stellar radii. 



2. Scaling to similar solutions 



Although the phenomenon of core collapse is quite complex (e.g., Colgate 1995| ), the resulting 



blastwave in the stellar envelope has relatively simple dynamics. Gravity is relatively insignificant 
except for the interior layers, which may suffer infall. The gain or loss of energy to nuclear reactions 
is also dynamically insignificant except for these inner layers, although they are probably mixed to 
higher velocities by the Rayleigh- Taylor instability. (Heat from decaying 56 Ni may blow bubbles in 
the ejecta [Arnett, Pryxell & Miiller 1989; Li, McCray & Sunyaev 1993]). Likewise, neutrino heating 
is negligible outside the dense circumnuclear material. The pressure is radiation dominated, and 
the ratio of radiation to gas pressure remains constant (for each mass shell) throughout the period 
of adiabatic motion. The radiation diffuses sufficiently little during the period of acceleration that 
the flow can be modeled as adiabatic, although this approximation fails for the optically thin shell 
associated with the UV or X-ray flash (e.g., Ensman & Burrows 1992). The heat lost to radiation 
never exceeds a few percent of the total explosion energy, implying that radiative transfer can be 
neglected in general. 



Under the approximation that the flow is adiabatic, radiation-dominated and non-gravitating, 
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its only dimensional scales are the explosion energy (E{ n ), ejected mass (M e j) and initial stellar 
radius (i?*). Explosions with different values for these parameters, but the same initial density pro- 
file, are identical after rescaling the hydrodynamic variables by their characteristic values. Common 
stellar density profiles recur on a variety of scales, and the resulting (normalized) ejecta distributions 
are identical. We will make use of this principle in §|5[ 

We take a mass coordinate m{r) = M(r) — M rem (where M rem is the mass of the material 
that will become the remnant) that is zero at the boundary between ejected and remnant material, 
and increases to M e j at the exterior. The configuration of the progenitor is specified by the initial 
density po(m) and the initial radius ro(m). If we normalize to the characteristic scales of mass, 
radius, velocity, time, density, pressure and entropy, the flow expressed in normalized variables 
evolves independently of these scales. The basic characteristic scales are Af e j, and E- m , from 
which it is useful to derive the following characteristic quantities: 

E- \ 1 / 2 ( M ■ \ - 1 / 2 



„ = _| _ 4.720 X 10- (^jyy gem-, 

» s fr^ 2 - 373 xioI3 (io^i)(^)"^^' 



We will call the quantity s = p 1 / 7 / p the 'entropy'; for 7 = 4/3, it is proportional to the phase space 
volume and the photon-baryon ratio. 

The goal of this paper is to provide formulae for the end product of an adiabatic explosion. 
In the final state of adiabatic free expansion, the velocity of a given mass shell is fixed, while its 
density and pressure fall off as t~ 3 and t -37 = i -4 , respectively. The ejecta properties scale with the 
ejected mass, stellar radius and blastwave energy: for a given feature in the ejecta (the half-mass 
shell, say), 

v f oc v* = Ell 2 M^' 2 

Pf oc P ,(t/u)- 3 = Er^ M yr 3 

p f oc p*(t/t*)- 4 = Er^M^r 4 (2) 

These overall scalings are very important for inferring gross quantities, but they say nothing about 
the relation between the progenitor's structure and the distribution of its ejecta. 

So, we are interested here in the normalized ratios Vf/v+, p/t 3 / ' p*t 3 , and pjt 4 /p+t 4 . as parame- 
terized, for instance, by the normalized mass m/M e j. Note that these are 'final' values only in the 



- 9 - 



sense of an adiabatic explosion into vacuum: the radiation pressure will differ from pft 4 because of 
radiation diffusion and radioactive decay, and the distribution of matter in the ejecta (vf(m) and 
Pf(m)t 3 ) will ultimately be destroyed as interstellar material is swept up. 

Throughout this paper, we use the normalized mass coordinate rh = m/M e j, which varies from 
to 1 in the ejecta. Because all of the hydrodynamical variables approach power laws in the limit 
of high velocity, rh — > 1, we introduce a notation that removes this dependence. In this limit, a 
variable x varies as x oc (1 — fh) Wx for some exponent w x ; so, we define f x (fh) by, 

x = xj x (m){l - m) w * , (3) 

where is a characteristic value for x (a combination of E- m , M e j and R±). The exponents w x 
and limiting coefficients f x (l) can be evaluated by combining the initial density and shock velocity 
in the subsurface layers with results from the planar, self-similar solution for the ejection of this 



material; this is done in §5.1. Unless otherwise noted, f x refers to the final value of x in free 



expansion. So, /„ refers to Vf, whereas f Vs refers to v s , the shock velocity. 



3. Stellar structures of the progenitors: stellar models and approximations 

The progenitors of core collapse supernovae are stars initially more massive than about 8M Q , 
whose cores have advanced though carbon burning. The degenerate iron core approaches the 
Chandrasekhar mass in the subsequent burning stages, the process of 'core convergence' (Arnett 
1997). This core has the radius and density of an iron white dwarf. The circumnuclear shells 
decrease in density away from this core, characteristically as po °^ r cT 3 - Outside the carbon shell, 
the helium envelope begins. There is typically a steep, shelf-like decline in the density between the 



carbon and helium. The helium envelope is typically radiative prior to the explosion ( Arnett 1997 ) 



Above the helium envelope sits whatever hydrogen-rich envelope has survived the stellar mass 
loss, which may be none at all. The structure of the star interior to the hydrogen is little altered by 
its presence or absence, because the hydrogen burning shell has a pressure and temperature that is 



negligible relative to the interior (Arnett 1997). In stars that have hydrogen, the boundary between 
hydrogen and helium is marked by a density drop even steeper and greater than at the base of 
the helium envelope. Presupernova stars with radiative hydrogen envelopes are blue supergiants 
(BSGs); those with convective hydrogen envelopes are red supergiants (RSGs); and those with no 
hydrogen at all are helium or carbon stars. 

Core-collapse supernovae are classified observationally into types II, lb, and Ic. Type II su- 
pernovae show hydrogen lines, type lb show helium but no hydrogen, and type Ic show neither. 
Because the spectrum may change with time, supernovae can shift in classification (Filippenko 
1997): supernova 1993 J was one such (e.g., Shigeyama et al. 1994). Thus red and blue supergiants 
have the capacity to produce types II and lb, helium stars can produce types lb and Ic, and carbon 
stars can only produce type Ic. 
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We will concentrate on the topmost stellar envelope, which we will term the outer envelope 
and whose mass we designate M env . This layer is the hydrogen-rich envelope of those stars that 
have hydrogen. Those layers that lie interior to the outer envelope, but are also ejected in the 
explosion, we will refer to as the mantle, with mass M man . The total ejecta mass M e j is the sum 
of these, M e j = M man + M env . The remnant, with mass M rem , is the stellar material that is not 
ejected; this is slightly greater than the mass of the relict object because of rest energy radiated in 
the explosion. The total progenitor mass is M* = M rem + M e j. Our focus on the outer envelope 
is motivated by the observation that it receives the lion's share of the explosion energy; moreover, 
our simulations indicate that the distribution of the outer envelope ejecta is determined primarily 
by its own structure and the mass M man of the ejecta it overlies, not appreciably by the structure 
of the underlying mantle material. 

Convective envelopes, like those in RSGs, can be approximated as n = 3/2 polytropes in 
structure. This approximation would be exact if the convection were efficient, but the presence of 
a superadiabatic gradient breaks this idealization at large initial radii (Nomoto & Sugimoto 1972). 
Models of convection can prescribe a superadiabatic gradient, but the results depend on the model 
used (mixing length, for instance), as well as stellar parameters (such as L^/L^). We will idealize 
convective outer envelopes as n = 3/2 polytropes, but the reader is advised to keep this caveat in 
mind. 

Radiative envelopes of constant (Thomson-dominated) opacity, such as those in BSGs, limit 
to an effective polytropic index of n = 3 toward the outside of the star, because such envelopes 
behave like n = 3 polytropes when luminosity and mass are proportional. However, the interior of 
blue supergiants is better described by a lower effective polytropic index. In the BSG progenitor of 
Nomoto, we find that 2.1 < n e g- < 2.4 for the inner 75% of the ejected mass, with n e s climbing to 
3 in the exterior 25%. (A power law of index 17/7 = 2.43 was adopted by Chevalier &: Soker 1989[ ) 



The structure of a polytropic envelope is determined, up to overall scaling factors, by a sin- 
gle parameter, corresponding to the fact that the Lane-Emden equation reduces to a first-order 
differential equation in terms of homology- independent variables ( Chandr asckhar 1 939 ) . Which 



structure is appropriate is determined by the inner boundary condition, defined by the mantle's 
mass ratio (M man + M rem )/M + and the radius ratio i? man /.R*. In order to get a single parameter 
from these two, consider that the Lane-Emden equation describing the outer envelope can be in- 
tegrated all the way to the center of the star. This hypothetical integration through the region 
actually occupied by the mantle will yield a polytrope that behaves as if it lay above a point mass 
at the origin. The value of the point mass compared to M*, which we term q, provides a single 
parameter with which to distinguish outer envelope structures. 

To understand the relation between q and the more familiar parameters (M man + M tera ) / 
and -R man /i?*, we must examine the central behavior of polytropes. Polytropes with n < 3 limit as 
r — » to the central distribution ( Chandrasckhar 1939| ). 



po oc {BRjr - If, (4) 
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for some constant B > 1 (equality when q = 1). In red supergiants (n = 3/2), one can typically 
ignore the envelope mass taken up by the volume of the mantle, both because the mantle radius 
is very small (-R man ~ and because the envelope mass is not highly concentrated towards 

the center {fa oc r ). So for RSGs, we can use the core-envelope mass ratio to parameterize the 
envelope structure: 

q ~ 1 - M env /M*. (5) 

In BSGs, an effective value of q is less easily determined, because of the variation in the effective 
polytropic index n e g . 

Polytropes with q = extend continuously and smoothly to the origin, where dp/dr = 
(because dp/dr = at the center). For any value of q greater than zero, the central density profile 
is p oc r~ n , as long as n < 3 (|Chandrasekhar 1939|) . In polytropic envelopes that are much more 
massive than their cores, so that q <C 1, this behavior is limited to a small region around the 
center (of similar mass), outside of which the structure returns to the flatter, q = form. On the 
other hand, envelopes with q <^ 1 are dominated by the gravity of their cores: the density profile is 
everywhere at least as steep as the central behavior. 



3.1. Outer density distributions of stars 

We will be especially interested in the region of a progenitor nearest the surface, because it is 
in this layer that the shock accelerates and casts away the highest-velocity ejecta. Fortunately, the 
equations of polytropic stellar structure are simplified for this material. The hydrostatic equation 
for a polytropic envelope can be integrated to give: 

i/n_ l/nfR* \ Q-^rnir) + M Tcra ]dr-^ 



where 



l n Lin i i \ -it, / c \ 

Po =Pi — - 1 Z=i ' W 



4% M*dr~ 



l/n _ GM* 

Pl - (n + l)iOV (?j 
and Po = Kpo- I n a shallow (outer) layer whose mass is negligible compared to M+, the ratio of 
mass integrals is approximately unity. Therefore, this region has an analytical density structure: 

Po = Pi(— -lT. (8) 



ro 

This is valid only for masses such that the enclosed mass can be approximated by M+; however, 
the region of validity may be relatively large in radius. For light envelopes whose entire mass is 
negligible (q = 1), this distribution holds for the entire envelope. Here, n is 3/2 for efficiently 
convective envelopes (although inefficient convection may invalidate this form) and exactly 3 for 
radiative envelopes of constant opacity (an analytical result in the outermost layers, where the ratio 
L/M is constant). 
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The stellar density coefficient p\ is related to the mass coefficient / Po (l), denned in equation 
), by the relation 



the corresponding exponent for po as a power law in (1 — m) is, 



*» = iTTT < 10 > 

For radiative envelopes of constant opacity, 

_ a(pm H ) i G 3 M3 (3 4 
Pl 1924 1-/3' 1 j 

where 1 — (3 = L^/L^d is the ratio of the stellar luminosity to the Eddington limit, or equivalently, 
1 — (5 is the ratio of radiation pressure to total pressure in the outermost layers of the star. From 
equation (^), for radiative stars, 

/M(1) .„. 148 ^v /4 ^mfj^v /2 , (12) 



M e J 1 / 4 \0.62J V1OM 



and for Thomson opacity, 



For convective envelopes with < q < 0.9 we find that, 

/ \ 2/5 

/po(1) - f (0.37 - 0.18^ + 0.096g 2 ) . (14) 

This fit to the Lane-Emden solution is accurate to within 1% for the given range of q. 

In summary, each layer of a star has a density distribution that is po °c r o n at its interior, and 
that may flatten out towards constant density before dropping precipitously near the edge of the 
shell. The appropriate value of n is 3/2 for efficiently convective envelopes and 2.1 — 3 for radiative 
envelopes. Besides n, the structure of a polytropic envelope depends on a central concentration 
parameter q; for the outer envelopes of RSGs, q is essentially the mass interior to the envelope, 
divided by the progenitor mass. There are sharp drops in density between the carbon and helium 
and between the helium and hydrogen shells; the latter is most drastic in the case of red supergiants. 



4. Model for shock propagation 



It is through this layered sequence of envelopes that the supernova shock propagates. A model 
for the shock velocity must account for the general deceleration as mass is swept up, punctuated 
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by accelerations as the shock encounters steep density gradients. The velocity of the blastwave 
shock is clearly affected by all of the material it has encountered, as well as the conditions at the 
shock front. But, the history of analytical shock velocity approximations has primarily been the 
search for a combination of local variables that controls (or at least mimics) the motion of the front. 
Formulae in the literature have typically given v s in terms of the local initial density po and radius 
ro, as well as the blastwave energy, E- m (an integral quantity). 



Kompaneets (1960) considered shock acceleration in the case of an initially exponential density 
distribution. He approximated the postshock pressure as a constant fraction of the mean energy 
density, leading in the case of spherical symmetry and in our notation to the expression, 



where Tk60 is a constant that can be evaluated in specific cases. As discussed by Koo & McKee 
(1990), this approximation overestimates the shock velocity as the shock accelerates, both because 
the mean pressure drops in relation to the mean energy density, and because the postshock pressure 
drops in relation to the mean pressure. Laumbach & Probstein (1969) developed an approximation 
for the motion of the shock wave in which the shocked material is assumed to be concentrated in a 
thin shell behind the shock, moving at the postshock velocity. The approximation involved Taylor 
expansions for the form of this shell, and related these to the motion of the shock. As Koo & 
McKee (1990) discuss, this approximation tends to underestimate the shock velocity, and fails for 
strongly accelerating shocks. 

Klimishin & Gnatyk (1982) considered shocks that alternate between periods of deceleration 
and acceleration, as shocks in supernovae do. They proposed that the shock velocity should obey 
v s oc e\1 2 {pqTq)^ 1 / 2 during periods of deceleration, and v s oc E^ 2 (poro)~^ 1 (with j3\ = 0.2) during 
periods of acceleration. 

Koo & McKee (1990), considering density distributions that steepen monotonically, proposed 
a shock velocity that accelerates beyond a characteristic radius. Specifically, they chose a postshock 
pressure that changes from the form appropriate to spherical, decelerating blastwaves (as in Sedov 
1959) and that appropriate to strongly accelerating blastwaves (as in Whitham 1974) at a transition 
radius R s ^- This leads to the shock velocity formula, 



4.1. 



Previous models 




(15) 



Vs r (pO 
— =1 KM90 — 



) 



-Pi 



) 



3/2 




) 



-fa 



S.t 



(16) 



ro + -Rs,t 



where TKM90 is a coefficient chosen to match Sedov- Taylor blastwave if the central density distribu- 
tion is flat. This prescription worked well for the relatively smooth distributions considered, which 
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were flat at the origin and fell to zero density asymptotically at large radii with a well-defined scale 
length. 

However, supernova progenitors have layers of different compositions, separated by shelf-like 
drops in density at the burning shells (§|3|). Therefore one must employ a model for shock accelera- 
tion that does not rely on the choice of a single radial scale. The ability to match the acceleration of 
the shock in the interior shell boundaries and at the outer edge of the star is particularly essential. 



4.2. Improved model for shock velocity 

During a period of self-similar shock propagation (e.g., Sedov 1959), the shock velocity obeys 
v s oc y/Ei n /m(ro). During this period, as well, the ratio p(ro)/p(ro) is a constant, where p( r o) is the 
average density of the material to be ejected within ro (note that ~p and fa do not include the stellar 
core mass that will become the remnant object). During a period of strong planar acceleration, 
when ro and m(ro) are essentially constant (and so 7i(ro) is too), the shock velocity takes the form 
v s oc [p(?"o)/p( r o)] _ ^ 1 • We propose a continuous and simple form. To accommodate both spherical 
deceleration and planar acceleration, we take the product form v s oc ^/i?i n /m(ro)[po( r o)/7'o( r o)] _ ^ 1: 

a Way* fa')--'*. (it) 



\p*J \R* 

Using the approximation that a strongly accelerating shock runs along a forward characteristic of 
the postshock flow, Whitham obtained ft ~ 1/(2 + v^tAt = 1)], or 0.2071 for 7 = 4/3. The 
self-similar, planar solution of Gandel'man & Frank-Kamenetsky (1956) and Sakurai (1960) gives 
ft = 0.1909 and ft = 0.1858 for n = 3/2 and n = 3, respectively, so ft ~ 0.19 is a reasonable 



approximation. (Note that in writing equations [15]-|l7j for v s /v* we have suppressed the explicit 



1/2 

dependence v s oc E^ that is important for blastwaves whose energy is not conserved.) 

It is possible to define T(tq) in terms of moments of the distribution po(r < ro) in such a way 
as to make eq. (|l7|) strictly correct for self-similar blastwaves. For a self-similar blastwave in the 
distribution po oc r kp , equation fll~7| ) is correct if V = cr _1 / 2 [47r/(3 — k p )]~P l , where a = E- m / {mv^). 
For the Primakoff blastwave, a = 8/(3(7 + l) 2 ]) so that, 

' ' (7n)~^ (18) 



x 24 

for 7 = 4/3. Radiative stars roughly resemble the Primakoff power law po oc r (Chevalier & 
Soker 1989), as do the mantles of RSGs. Despite the actual differences between stellar distributions 
and this power law form, we find that together, equations (|17]) and ( |i~8| ) are remarkably accurate 
(Figure |). 

Throughout the bulk of the stellar mantle and envelope in our simulations, the exact choice 
of ft between the self-similar or approximate values makes very little difference. However, the 
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exact values are important for determining the power laws of high- velocity ejecta (§|5.1|). We find 
typical errors of v s in the stellar envelope to be ~ 2%, although higher (~ 10%) in the mantle and 
in regions of strong acceleration. But, these errors are in comparison to simulations of adiabatic, 
non-gravitating point explosions in a gas sphere of mass m{R+) = M* — M rem . In real supernovae, 
the shock velocity will differ significantly from this formula at the base of the ejecta, due in part to 
the mass that ultimately enters the remnant; in particular, the shock velocity does not approach 
infinity as m — > 0, contrary to equation (|l~7|) . Our focus is on the envelope and outer mantle, 
however, and for these equation (|l7]) should remain accurate. 

Taking the value 0i ~ 0.19 to represent both radiative and convective envelopes, we can write 



equation ( |17|) as: 



, , x -0.19 / _ \ -0.57 

v< ~ 0.794^m- a31 ^) ^ 



p*J \R* 
Wo/ 



0.794 ^) (^L) , (19) 
where m is the enclosed ejecta mass. 



To summarize the shock velocity approximations: the Kompaneets approximation has v s oc 
E 1 ! 2 p§ 1 ^ 2 r 3 ^ 2 , corresponding to a constant ratio of postshock to mean pressure. The Koo & 
McKee (1990) approximation goes from v s oc E^ 2 Pq^Tq 3 ^ 2 to v s oc E^ 2 p$ 131 r^ 2 at a transition 
radius i? Sit . Klimishin & Gnatyk took v s oc E^ 2 (po r o) -1 ^ 2 an( ^ v s ^ -E i 1 n //2 (po r o) _/31 m periods of 
deceleration and acceleration, respectively. Our approximation has v s oc E 1 ^ 2 m /3l ~ 1 ^ 2 p r 3/3x . 
By introducing m and using po and vq only in the combination pqTq, we avoid having to choose 
dimensional scales. We compare the Kompaneets, Klimishin & Gnatyk and Koo & McKee approx- 
imations to our equation (|i~7|) in Figure (S). 



Equation ([I?]) represents an improved approximation for v s in terms of the initial quantities m, 
ro and po and E- m . These Lagrangian variables are especially convenient because they are readily 
available in the initial state. Including m allows the dimensionless ratio Po(^o)/Po( r o) to enter 
into the shock velocity approximation. Whereas the shock velocity must in reality depend on the 
detailed structure of the material encountered so far, it appears as if this single moment of the 
initial density distribution is sufficient to specify the shock velocity quite accurately. 

It is important to note that equation (|i~7| ) is a purely local theory for the shock velocity, in 
the sense that it only makes reference to the variables (m,ro,po) a t the same mass shell in the 
progenitor. 

In the upcoming sections, we will be interested in relating the ejecta pressure to its density by 
means of the 'entropy', s(m) = p 1 / 7 /p. Using eq. ( [T7D with 7 = 4/3, we can specify the entropy of 
material that experiences no reverse shock: 



s(m) _ 6 3 / 4 3/2 r,„. ( ,-„r -"+«"•" ' -,..(,;,, - {>i < 2 



~Pq{™)~ 


-(1+6/30/4 


> (m)" 






[ R* \ 



^(6/3i-3)/4_ (2Q) 
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This entropy formula, like the shock velocity formula from which it is derived, is a local theory, 
invoking only (m, ro,po) f° r the- same mass shell of the progenitor. It falls short for shells that 
will encounter a strong reverse shock, a weakness that will prevent our models from describing the 
mantle ejecta accurately. Including reverse shocks would require more information than the local 
variables (m,ro,po) provide. 



5. High-velocity ejecta 



The analytical outer density distribution described in §3.1 provides a stable platform from 
which to launch a study of envelope expulsion in supernovae. The fact that polytropic envelopes 
of the same index share this common form is useful for understanding how a planar, self-similar 
analog is applicable to the shock and rarefaction flow in this region. Moreover, this distribution 
allows a characterization of the way in which the real, spherical flow differs from the planar analog. 

For material that is close enough to the stellar surface, (i?* — ro <Cfl*), the density distribution 
of equation (|8|) takes the limit p$ — > p\(l — ro/i?*) n . If a mass shell stays in this small range of 
radii throughout its period of acceleration, spherical curvature will not affect its dynamics. Then, 
its velocity will be set by the planar, self-similar problem first posed by Gandel'man & Frank- 
Kamenetsky (1956) and solved by Sakurai (1960), which we reconsider in the Appendix. The 
material covered by this solution settles into power law distributions whose indices are easily derived 
once the shock velocity index f3\ is known (c.f. Chevalier &; Soker 1989| ). If the coefficients of density 



[pi or /p (l)] and shock velocity are known, the coefficients of the final density distributions can 
also be derived (c.f. Imshennik fc Nadezhin 1989| ). In §5J. the indices and coefficients of this 



high- velocity material are calculated, using the results of §|| to set the shock velocity coefficient. 

This planar, self-similar analog is limited in its applicability to real supernovae by several 
effects. For mass shells that differ appreciably in radius from R±, po is no longer a simple power 
law in depth. For such shells, changes in radius are also sure to make a planar approximation poor. 
In fact, as first mentioned by Litvinova & Nadezhin (1990), curvature is actually negligible only 
when (1 — ro/ii*) 1 / 3 <C 1, not just when (1 — ro/R+) <C 1, because acceleration is very slow in the 
planar solution. 

The postshock and rarefaction flow are profoundly affected by sphericity, but in a way that is 
consistent among different progenitors. The effects that corrupt the planar solution are primarily 
geometrical in nature - that is, they occur because of differences in the initial distribution (with 
scale length i2*), or differences in the hydrodynamics (also with scale length i?*). Insofar as the 
initial distributions of material are identical after rescaling (i.e., so long as as eq. § is correct), the 
fully spherical processes of shock propagation and expansion in different progenitors are identical. 
This implies that the outer ejecta distributions of polytropes share common forms, which limit to 
the self-similar power law solutions at the highest velocities. Specifically, one can multiply these 
power laws by a function that depends on the initial depth, in order to account for the effects of 
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progenitor structure and spherical expansion. The validity of this approach is demonstrated in §5.2 
where such correction factors are also presented. This approach is similar to the analytical solution 
of Kazhdan & Murzina (1992), who present a correction factor (parameterized by 7 and n) that 
accounts for the first-order effect of sphericity at all times in the flow. As shown in the Appendix, 
the correction factors we derive resemble the outcome of Litvinova & Nadezhin's (1990) suggestion 
that spherical effects truncate planar expansion once the radius has expanded by a certain amount 
(of order R*). 

To put it another way, it is possible to define a set of variables that remove the scalings 
of physical quantities in the outer edge of a star. In terms of these variables, polytropes of the 
same index assume the same structure (no longer up to a multiplicative constant) at large initial 
radii. And, since the new time and velocity variables have taken out variations in the strength 
of the shock as it enters the outermost regions, different polytropic progenitors play out the same 
behavior of shock acceleration and outer envelope ejection - a behavior that includes the effects 
of the finite radius on the initial conditions and on the hydrodynamics of spherical expansion. 
Equation (|8|) demonstrates that the normalized initial radius tq/R* (or the normalized initial depth, 
xq = (R* — ro)/i?*) and the rescaled mass coordinate (1 — m)/(pi/p*) are variables of the required 
type, as is the density variable pjp\. (The results of the Appendix, as well as the analytical 
development of Kazhdan & Murzina [1992], show that the variable x^ 3 is more relevant than xo, 
and so we define the corresponding mass variable 5m = [(1 — m) / \pi / p*)} 1 / 3 ( n+1 ).) For the shock 
strength, we must rely on the results of §H|to generate a new velocity variable (v/v i[ )(pi/p^)^ 1 that 
removes the scaling of shock velocity with the structure of the star. In terms of these variables, 
the final distribution of high-velocity ejecta from polytropes of the same index n share a single 
distribution, which need only be determined once. Since it is convenient to write this distribution 
as a power law (the self-similar result) times a correction factor, it is natural to express the correction 
factor in terms of x l J 3 , 5m, or (vf/v i ,)(pi/p i ,)^ 1 , as is done in § |5.2| . 

However, such correction factors must themselves fail to describe the ejecta at some depth, 
or for velocities that are too low. This follows from the fact that equation (|8|) is valid only in the 
region for which the enclosed mass is essentially M*. For radiative envelopes of constant opacity, 
the effective polytropic index deviates from 3 for the interior mass shells. For polytropes as well, the 
initial density distribution is affected: the actual density distribution (equation ||) differs from the 
common distribution (equation ||) when m(r) + M rem < M*. Even when the shock wave is in the 
outermost mass zones, its velocity is still being affected by the forward characteristics of the earlier 
blastwave flow. At what point does the flow forget its history and limit to a common distribution in 
these outer regions? Theory (Whitham 1960) suggests that a shock's velocity behavior is determined 



by local conditions once it begins to accelerate strongly. The approach adopted in §5.2 is to 
determine numerically what region of the ejecta is reproduced by a variety of progenitors, and to 
present fits for correction factors only in this region. 

There is an upper velocity limit on the validity of our solutions as well, because the finite 
width of the radiative shock becomes comparable to the depth at some point. Our assumption of 
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adiabatic flow must break down by this point. We estimate the maximum velocity in §5.3 



5.1. Limiting power laws for the highest-velocity ejecta 

From our shock velocity approximation and the resulting entropy distribution, equations ( |i~7| ) 
and (|20|), the outer coefficients of shock velocity and entropy are, for 7 = 4/3, 

f Vs (i) = r/ po (i)-^, 

/ s (l) = 0.1273r 3 / 2 / w (ir (3ft/2+1/4) - (21) 

In the planar, self-similar solution, Vf{m)/v s {fh) is a constant, (vf/v s ) p , whose value we derive in 
the Appendix and present in Table [2]. This ratio gives f v (l) from f Va (l); then, mass continuity and 
entropy conservation give / p (l) and / p (l): 



fv(l) = /«.(!) (^) , 

\ v s / p 



-3 



s / p 

f p (l) = 0.002191 (^-) 4/3r_2 /po(l) 2/?1 - 1/3 (^)~ 4 - (22) 

The exponents of these variables in (1— m) are easily derived, but we state them here for convenience 
(the values are given in Table Q): 

(77/ 
--. 
n + l, 

(1 \ 2 n 

w s = - - 1 )w Po + -w v = (1 - 2/?i - 7), 

\7 / 7 7( n + !) 

77 



77J p = 1 — 3777,; = 1 + 3/9l 
777 p = 7(777 s + W p ) 



n + l, 
7 + n — 2p\n + 3/9i7n 
n + l 



(23) 



As an example, the asymptotic velocity distribution is: 
W/ (m-»l) = f v (l){l-m) w "v* 

= t(^) / P0 (l)-^l-m)^ in/(n+1 V (24) 

\ V S / p 

The asymptotic final density and pressure distributions are given in terms of velocity by pj oc 
v l j 2 t~~ 3 and oc vJ 2 t~ A , where 

, , 71 + 1 + 3/3in 

tp2 = Wp Wv = 7, , 

pin 
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7 + n + /?m(3 7 - 2) 

<p2 = w p w v = . (25) 

Pin 

(The subscript 2 denotes power laws for the outer ejecta - see §0.) The values for n = 3/2 and n = 3 
are given in Tables ||| and || respectively. Based on comparison with actual stellar models (§|9|), 
there is some indication that even when the local value of n e g = d log po/d log xq varies, equation 
( p5| ) still roughly approximates dlogpf /dlogvj and dlogpf /dlogVf for the high- velocity ejecta 
of peripheral material. Note that a shallower initial density produces a steeper final density and 
pressure in the outermost ejecta. 



5.2. Common spherical dynamics of high-velocity ejection 

In Figures || and |4] we demonstrate that convective and radiative stellar envelopes indeed each 
assume a common exterior distribution, as discussed above; that the shock velocity approaches a 
common form for each (up to a constant); and that the ratio v /(?"o) /v s (tq) does as well. In the 
Appendix we investigate the implications of the suggestion of Litvinova & Nadezhin (1990) that 
the final velocity is set by the truncation of planar, self-similar acceleration once a shell's radius 
has expanded to a radius somewhat larger than R+. This, along with the form of the self-similar 
dynamics (equation |A5[| ) gives a form for vj/v s , equation (A7) that is a polynomial in x^/ 3 ', where 



xq = 1 — vq/R+ is the normalized depth. It is no surprise, then, that the numerically determined 

1 /3 

distribution of Vf(ro)/v s (ro) is best fit and expressed as a function of x (as shown analytically 
by Kazhdan & Murzina 1992): 

v _fM „ (H) (i _ o.SIxy 3 + 0.76x 2 J 3 - 1.19x0) (n = 3/2), 
v s (x ) \v s J p 

Yl^l ~ (Hi) (1 _ 0Mxl /3 + 0.24^ /3 - 0.37x ) (n = 3), (26) 
v s (x Q ) \v s J p 

where (vf/v s ) p is given in Table ||. These fits were made to the outer half of the radius of the 
q = 0.3 polytrope in each case. The fit for n = 3/2 is good to within 1% for tq > 0.55 for both 
q = and q = 0.3 and to within 3% for > 0.2 for q = 1 polytropic envelopes. The fit for n = 3 
is good to within 1% for vq > 0.5, 0.6, and 0.65 for q = 0.3, q = 0, and q = 1 respectively. As 
shown in Figure [| the full spherical result, equation (^), resembles the outcome of Litvinova & 



Nadezhin's suggestion, equation (A7), if acceleration halts at about 3-R*. 



We note that the first-order correction terms in equation |2^ are predicted by the analytical 
perturbation theory of Kazhdan & Murzina (1992). From their theory, the first-order coefficient 
should be —0.3253 for n = 3 above, where we find —0.34 in our fit. The two coefficients are 

somewhat different in meaning, in that the Kazhdan & Murzina coefficient represents a term in 

1/3 

the expansion of Vf/v s around x =0, whereas ours is part of a cubic polynomial fit to a range 
1/3 

of the variable x . The value —0.3252 for the first term in this expansion implies that planar 
acceleration of the outermost material is effectively truncated at 2.94i?* in Litvinova &: Nadezhin's 
scenario, very close to our estimate of 3i2* (see the Appendix). 
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In the region for which envelopes of the same value of n scale to a common solution, we can 
give the Lagrangian form of the ejecta distributions, i.e., Vf(m), pf(m)t 3 and pf(m)t 4 . We present 



these high- velocity formulae in Table 0. These formulae consist of the power laws derived in §5.1 



multiplied by a function of mass that accounts for the deviations from self-similar flow. Since the 
mass external to a depth xq is 1 — m(xo) = 4ttxq +1 (pi/ p±)/(n + 1), the mass variable corresponding 
to Xq^ 3 can be defined as: 

5 m ^[(l_^) (/0l /^)-l]l/3(n+l). ( 27) 

This variable is convenient for giving the correction factors to power laws in (1 — rh), as in Table 

rt 1 /3 

^, because the effects of spherical expansion are linear in 8m just as they are in x ' . 



5.3. Shock breakout and the maximum ejecta velocity 

When a supernova shock approaches the surface of the star, there comes a point at which the 
postshock radiation can leak out in a burst of ionizing radiation (e.g., Chevalier & Klein 1979). 
This event marks the end of shock acceleration, and limits the velocity at which ejecta can be cast 
away. The shock velocity at breakout is well-determined from simple considerations. This sets the 
upper velocity limit for the assumption of adiabatic flow, and the gross properties of the radiation 
outburst ( Imshennik &: Nadezhin 1989| ). 



Shock propagation is unaffected by the leakage of radiation until the optical depth of the shock 
exceeds the optical depth of the available material. Since the shock wave is a region in which light 
diffuses upstream and gets amplified by compression, the time for light to diffuse across the shock 
front must be comparable to the time for material to flow through it. This equality sets the physical 
size of the shock transition: a simple estimate shows that the shock's optical depth is 

r s ~ -. (28) 

v s 

This is only an approximate equality, because the shock structure set by radiation diffusion is 
not sharp. In fact, the structure of a strong, radiation-dominated shock front propagating into a 
uniform medium can be calculated analytically using the diffusion approximation (Weaver 1976). 
In this solution, the radiation pressure varies with optical depth as: 

P = 6exp [-3(t -t )v s /c] 
p v 2 s 1 + 7exp [-3(r - t )v s /c] ' 

where t is the optical depth (increasing inward) and To is a constant that defines the position of 
the shock front. Here, the radiation pressure approaches its upstream and downstream asymptotes 
exponentially. However, the smoothness of the shock front introduces little uncertainty into the 
velocities (~ 20%) and temperatures (~ 10%) derived below ( Imshennik fc Nadezhin 1989| ). 



Shock breakout occurs when r s is equal to the stellar optical depth t(xq) at some depth xq: 

r(x ) = Kpl *f[ +1 , (30) 
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or, equivalently, when the radiation diffusion time across the shock front equals the time for 
the shock to get to the stellar surface (the local dynamical time). The terminal shock velocity 
t> s (breakout) at this point, where r s = t(x ), is, from equations (|i~7|), (28), (p0|): 

(31) 



v s (breakout) ~ 



-l/n 



'kM p 



1 + nJ \c ) V Rl 

The breakout shock velocity scales most significantly with v±, because (3\ is a small number. 

It is interesting to note that the width of the shock at the time of breakout is about 0.014i?* 

1/3 

and 0.019i?* for RSGs and BSGs, respectively. So, the parameter x is at least 0.24 and 0.27, re- 
spectively, for the adiabatic ejecta. Therefore the planar, self-similar limit is never strictly achieved. 



5.3.1. The speed limit for supernova ejecta 

One might expect that after breakout, the diffusion of radiation would interfere with adiabatic 
expansion at lower and lower velocities. If this were true, even those mass shells not involved in the 
shock at the time of breakout would be affected, and our adiabatic results would not be applicable 
to these shells. 

However, a straightforward comparison of the diffusion time (\^3t 2 /cup) and the dynamical 
time [t — t(breakout)] for a particular mass shell shows that this does not actually happen. Instead, 
the adiabatic approximation holds for mass shells deeper than a few shock widths from the surface 
at the time of breakout. 

To see this, consider the phase of planar free expansion, when a mass shell has traveled far 
compared to its initial depth, but only a small distance relative to R*. Then, because the surface 
area of a shell is effectively constant, the density scales as p(m) cx l/[t — i(breakout)]. Since the 
optical depth r of a mass shell does not vary in planar flow, the diffusion time is proportional to 
the dynamical time during this planar free expansion phase. Examining the planar, self-similar 
solution presented in the Appendix, we find that the ratio of the local diffusion time to dynamical 
time is never less than its value in free expansion. Similarly, the shock front itself has a width that 
is set by the equality of the diffusion time across the front to the time for material to cross the 
shock front. As a result, the criterion for a mass shell to experience significant radiation diffusion 
during planar expansion is simply that it be entrained in the shock at the time of breakout. 

Even so, one might still worry that the onset of spherical expansion will permit a diffusion 
wave to propagate inward. Then, the geometric reduction of a mass shell's velocity would be 
enhanced (and the formulae of § |5.2| invalidated) by the loss of photons. Fortunately, this does not 
happen either. The results of § |5.2| indicate that the geometric effect sets in at the fixed radius 3R+. 
Expansion to a radius of 3i?* decreases the ratio of diffusion time to dynamical time by only a 
factor of about ~ 9 relative to its planar value. One only has to move to slightly (~ 30%) lower 
velocity to recover this factor between the two time scales. Therefore, geometrical expansion ends 
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adiabatic acceleration before radiation diffusion ever becomes important - except for the material 
caught in or very near the shock front at the moment of emergence. (See Chevalier 1992 for an 
analysis of a diffusion wave in spherical free expansion.) 

So, we are justified in applying the full, adiabatic, planar acceleration factor (vf/v s ) p to the 
shock velocity at breakout in order to estimate the upper velocity limit Vf ma , x for which adiabatic 
results hold. 

For convective (n=3/2) and radiative (n=3) envelopes of typical RSG and BSG dimensions, 
respectively, we find: 
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Note the insensitivity of Vf max to the structure parameter pi/p*. 



(32) 



5.3.2. Relativistic mass ejection 

Recently, Kulkarni et al. (1998) have inferred that the radio-bright shell around the peculiar 
supernova 1998bw expanded relativistically during the period of observation, and take this to be 
physical evidence of the possible association (Soffita et al. 1998) with the gamma-ray burst GRB 
9804025. Woosley, Eastman & Schmidt (1998) considered models for 1998bw with ejected masses 
in the range 5-12 M & and energies in the range 4-30 xlO 51 ergs. We are motivated by these 
developments to reconsider the question of relativistic mass ejection, addressed initially by Colgate 
& Johnson (1960) and Colgate &: White (1966) and most recently by Woosley et al. (1998). 

Our criterion for relativistic mass ejection is that the maximum ejecta velocity given by our 
non-relativistic theory, equation (^2[), exceed the speed of light. (For the behavior of shock waves 
in the relativistic regime, see Colgate & Johnson 1960 and [Johnson fc McKee 1971| ). The star must 



be sufficiently compact to meet this criterion. For instance, applying equation (^) to an explosion 
of 10 51 ergs in a white dwarf, t>/ max would be of order 1.7c if n = 3/2, climbing to ~ 5.4c if n = 3.25 
in the atmosphere. But, the helium star progenitors of ordinary type lb and Ic supernovae may 
not be compact enough. Taking stellar and ejected masses for such progenitors as considered by 
Woosley, Langer & Weaver (1995), and using their radii at carbon ignition (1-9 Rq) and an effective 
outer polytropic index n e g = 3 to represent the presupernova star, explosions of 10 51 ergs produce 
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maximum final velocities in the range 0.5c-0.7c. 

In the periphery of a star, the shock velocity is a simple power law of initial density (equation 
|jl7f ), with the coefficient given by equation (|l8|). In order for the non-relativistic extrapolation 
of the final velocity of a mass shell to exceed c, its initial density must be smaller than the mean 
density by a critical factor: taking f3\ = 0.19 and Vf/v s = 2.1 in equation ([17|), 

Po & P*. (33) 



But, the electron-scattering optical depth of this material must exceed about 2, or the radiation 
driving the shock will escape. For an outer density law p$ = Pi(l — r o/R*) n , this criterion puts an 
upper bound on the radius: 
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If the star is sufficiently compact by the above criterion, equations (|9|) and give a rough estimate 
for the mass of relativistic ejecta: 

10 29.6-4.8/n Z^N-l/" / E - m \™(W/n) , M X 1-2.6(1+1/,) 

Mr ^^Ti— y ii^j It^j g - (35) 

We have assumed that the initial depth of this material is small compared to R*. An explosion 
of 10 52 ergs ejecting 1M & in a sufficiently compact star (with n = 3 or, as considered by Colgate 
&; Johnson 1960, n = 3.25) will produce relativistic ejecta with M re ic 2 ~ 3 x 10 48 ergs, for typical 
values of pi/p+. (But, to compare, we would predict a shock velocity 50% higher than Colgate & 
Johnson found, at the point where po = 30 gem -3 in their progenitor.) This energy is comparable 
to the minimum energy associated with relativistic electrons in SN1998bw (10 49 ergs; Kulkarni et 
al. 1998), and greater than the energy in the gamma-ray burst (8.5 x 10 47 ergs; Galama et al. 
1998), if the burst and supernova are associated. 



5.3.3. Observable aspects of shock emergence 

The shock velocity approximation of §|] allows us to predict the variation of the breakout 
quantities with the explosion parameters E- m , R±, M* and pi/p*. The quantities that can be 
extracted from simple analytic calculations include the postshock radiation temperature T se , the 
outburst energy .E se , and the time scale on which radiation diffuses out of the shock, t se (assuming 
non-relativistic ejection). Detailed information about shock emergence, like the hydrodynamics 



of optically thin material (which has been suggested to form a thin, dense shell, e.g., Ensman 



fc Burrows 1992 ) or the spectrum of the outburst, are beyond the scope of our simple analysis. 
The outburst spectrum is affected by the work that photons do as they diffuse out of the shock; 
a detailed solution awaits multifrequency radiation hydrodynamics simulations. Nonetheless, the 
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values of T se given below agree with recent estimates of the peak color temperature in outburst 
(e.g., Ensman fc Burrows 199^ ) within 10%. 



The postshock radiation temperature is related to the postshock pressure at the time of shock 
emergence, so aT^/3 = 2 / oo^s( Drea kout)/(7 -f l) : 
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The energy E se of the radiation outburst can be estimated as the thermal energy in the shock 
front at the time of breakout. So, E sc ~ (aT^/3)[47rit^xo (breakout)]: 
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The energy E se will be released on the diffusion time at shock breakout, which is also the time 
for the shock to travel its width. So, 
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However, note that the pulse will be longer than this from the vantage of a distant observer because 
of the light travel time (e.g., Ensman &: Burrows 1992] ), which is 1160 seconds for a star of radius 



500i? Q , or 116 seconds for a star of radius 5Oi? . As a result, the observed luminosity will be 
typically somewhat less than E se /t se . 
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These formulae reproduce the energy and time scale of the radiation outburst in the numerical 
simulations of Ensman & Burrows (1992), for explosion parameters chosen to match theirs. These 
formulae are also in agreement with the analytical results of Imshennik & Nadezhin (1989) for 
SN1987A, but only if we account for the fact that the shock velocity they use is 47% faster than is 



given by equation (17) for the same model, as if V were 1.16 instead of 0.79 for the outermost mass 
shells. If instead we use equations (^) and (j37[), we find a value of T se that is 20% lower, and a 
value of E se that is 50% lower, than those of Imshennik &: Nadezhin. It is not apparent why the 
two shock formulae do not agree, as their formula is taken from their simulations, and our formula 
agrees with our simulations. 

For the helium star progenitors of typical type lb and Ic supernovae, described by Woosley, 



Langer and Weaver 1995 and discussed in §5.3.2 (again, assuming that E in = 10 51 ergs), the energy 
associated with the outburst would be between 3 x 10 and 2 x 10 46 ergs, and the outburst would 
last between two and twenty seconds. In making these estimates, we have assumed that the stellar 
wind does not have sufficient optical depth to support a radiation-dominated shock. For compact 
progenitors with dense Wolf-Rayet winds, this assumption may fail. In that situation, the shock 



velocity formula (equation 17) can be applied to the stellar wind, and breakout quantities estimated 



at the radius for which the wind optical depth matches ~ c/v s . The circumstellar interaction (e.g., 
Fransson, Lundqvist & Chevalier 1996) would then begin immediately. 



6. Pressure-based model for the ejecta pressure and density distributions 

We now wish to present a model that describes the distribution of all the ejecta, not just the 
high velocity ejecta discussed in §g. To do so, we will extend the form of the pressure distribution 
inward in mass, approximating its variation with a simple functional form. In our calculations, the 
final pressure distribution is invariably the smoothest of the hydrodynamical variables (expressed 
as functions of rh) ; this is not surprising, since the pressure gradient is inhibited on small scales by 
its ability to accelerate the material so as to reduce its magnitude. Once we make a model for the 
final pressure distribution, the final density distribution follows immediately from the entropy left 
behind by the forward shock. 

To make a simple model for the final pressure distribution, py(m)i 4 , we multiply the pressure 
distribution known for the high- velocity material, given in Table (||), by a simple function of m. 
This has the advantage of preserving the high-velocity behavior of the flow derived in §|5], because fh 
varies very little in the region of validity of this solution. The solution found in this manner is less 
accurate than the solutions of §|5|, but they have the advantage that they are robust to variations 
in the progenitor structure (like the existence of superadiabatic gradients in RSGs). They fail in 
regions that experience a strong reverse shock, especially the mantles of RSGs. However, as we 



discuss in §6.1, the density jump between the mantle and outer envelope ejecta can be predicted 



despite the formation of reverse shocks. 
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Multiplying the pressure distribution specified in Table || by a simple function of mass produces 
our model for the final radiation pressure distribution: 



Pf(rh)t 



Pf(rh)t 
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4/3 
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(39) 



where [p/(w)t 4 /p^tf] pL+HV is the pressure distribution given in Table [| that is, a self-similar outer 
power law times a correction factor that accounts for the geometrical effects spherical expansion. 
The new term, [1 — a(l — m)] 4 / 3 , is but another correction factor that accounts for our imperfect 
knowledge of the pressure distribution in the bulk of the ejecta; we have chosen the power 4/3 to 
simplify the result for the final density distribution. It only remains to specify the free parameter a 
used in formulae p9|), fl4l"D and fl42]). In §|8]we will develop a technique for generating realistic RSG 
progenitors with which to constrain a; the variation of a with the stellar structure parameters of 
red supergiants is presented in Table ||[ 



The final density distribution corresponding to the ejecta pressure model of equation (|3 
obtained by means of the entropy deposited in the forward shock: 
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where s(m) is given by equation (20). Writing out the resulting density model in terms of the 
variables m, tq and po, we find 
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The only parameter remaining to be determined is the numerical prefactor A p . The planar, 
self-similar solution gives the definite value A p = ll/ p (l) 3 / 4 , where f p (l) is given by equation (f22|). 
However, we find that using overall energy conservation to set A p gives more accurate results. To 
do this, we first derive the velocity distribution from the density model by means of a numerical 
integration: 
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(Note that by integrating from zero velocity, we have excluded the remote possibility of hollow 
blastwaves, which can occur in special distributions; see Ostriker & McKee 1988.) If we define the 
model energy ^modei by, 

ffmodd = 1 f 1 v f (m) 2 
E- m 2 



-dm', 



(43) 
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then the energy of the ejecta model can be brought into agreement with E- m by the transformation 
(from eq. [[]): 
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In this transformation, the characteristic scales (v+, etc.) are kept fixed but the model ejecta 
distributions (vf(rh), etc.) are altered to agree with the blastwave energy. 

Equation ( [4l|) constitutes a local model for the final ejecta distribution: one that specifies the 
final density of a mass shell by reference to its initial variables (m, ro,po) only. The enforcement 
of overall energy conservation, and the application of the parameter a, weaken this statement only 
slightly. This model inherits from the entropy and shock-velocity models its inability to describe 
regions that have been hit by a strong reverse shock. An underestimate of the entropy leads to an 
overestimate of the density in the mantles of supergiants, RSGs especially. 



6.1. Density variation across the outer composition boundary 

Independently of our model for the distribution of the outer envelope ejecta, our shock velocity 
model allows us to evaluate the relative jump in density at the base of the outer envelope. This 
boundary is typically the outermost composition boundary in the progenitor. If we take pf{m)t A 
to be continuous in rh across the boundary, the entropy left behind by the shock wave, equation 



d4l|) with /?i = 0.19, dictates: 

5log( Pf t 3 ) = 0.545 log(po) +0.865 log(ro), (45) 

where 5log(po) and <51og(r ) are to be evaluated for nearby mass shells straddling the composition 
boundary. 

One might be concerned that this result would be invalidated by the formation of a reverse 
shock at this boundary, or by the development of the Rayleigh- Taylor instability. However, in our 
simulations the compression wave that decelerates the mantle material travels inward (in mass) 
before it steepens into a shock wave, so that the entropy distribution given by equation ( ^0[ ) 
is correct across the boundary. The Rayleigh- Taylor instability may indeed break the spherical 
symmetry of the ejecta. However, it is unlikely to change the fact that the outermost mantle 
ejecta is in contact with the innermost ejecta of the outer envelope. Then, pressure equilibrium 
across the boundary gives the density jump above. An interesting corollary is that the temperature 
distribution is continuous as well, even once radiation diffusion becomes significant. This may 
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have implications for the recombination epoch, since helium (with a much lower specific entropy) 
underlies any hydrogen-rich ejecta. 

Those interested in modeling the distribution and expansion velocity of the radioactive inner 
mantle material should note that our model for the outer envelope ejecta constrains these quantities 
rather poorly, because we overestimate the density of the mantle material. 



7. Harmonic mean models for ejecta distributions 



Our simulations demonstrate that the final density and pressure in the hydrogen ejecta can 
often be fit reasonably by the simple density model, 
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These models interpolate between limiting power laws at low and high velocities, joined at a specified 
value (/Obrcak or Pbrcak) at a given break velocity (f p break or v p break)) with a curvature set by the 
parameters y p or y p . 

For supernovae in red supergiants, these models apply only to the hydrogen ejecta, because of 
the severe density jump at the base of this material. As a result, we cannot apply the constraints 
of mass and energy conservation to give relations between the five parameters p&i 3 , v p breaks Ipi, lp2 
and y p in the density model. Instead, we fix l p 2 and l p 2 at their self-similar values; we set y p and 
y p at characteristic values, and we fit the other three as functions of q, using the techniques of §|8|. 
Table || presents these fits for red supergiants. 

These simple models fit the outer envelope ejecta of polytropic models for red supergiants quite 
accurately, to about 9 — 11% error in density and 23 — 30% error in pressure (both errors increasing 
with increasing q). In fact, the residuals have a very consistent pattern in each case, indicating 
that more complicated models could achieve even greater accuracy. However, real stars may not 
resemble polytropes sufficiently to justify any more complexity. 

For supernovae in stars with radiative outer envelopes, such as blue supergiants, no sequence 
of simple models models is available for the fitting of these parameters. However, the difference 
between mantle and outer envelope is not as severe for these stars, so we can apply models to the 
entire ejecta at once. The constraints of mass and energy conservation then reduce the number 
of parameters in the density model from five to three. The outer power law behavior is strongly 
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constrained by the planar, self-similar solution. Chevalier & Soker (1989) conjectured that the 
interior power law could be described by self-similar rarefaction, a hypothesis we also adopt to 
set I p i and l p \. Taking y p to be the same as in the RSG models, the density distribution is fully 
specified. For the parameters y p , v p break and Pbreak of the model pressure, few constraints are 
available. We quote typical values only. The result is a fiducial model, given in Table ||, that may 
be representative, but which depends not at all on the actual progenitor structure. In this sense, 
the pressure-based model of §|6| is superior, especially for radiative progenitors. See Table ^ for a 
synopsis of the relative accuracy of these models for realistic blue and red supergiant progenitors. 

8. Survey of ejecta distributions from model red supergiants 

The extreme difference in density between the mantles and outer envelopes of red supergiants, 
as well as the relative constancy of the hydrogen envelope density, suggest that simple models 
can reproduce the dynamics and final density distribution of these stars. Such an approach was 
taken by Chevalier (1976), who found that progenitor models with constant-density mantles and 
constant-density outer envelopes, separated by a region of steep density decline, could reproduce 
the qualitative features of type II light curves. Most notably, he identified existence of a plateau 
phase with the formation of a 'detached shell' (density maximum) in the final distribution. We will 
use a different set of simplified models to represent these progenitors, so it is worthwhile to review 
in detail why such models are realistic. 

In red supergiants, the outer convective envelope can be modeled as a polytrope of index 
n = 3/2 (§||), so that the shape of its density distribution is determined by the ratio of its mass to 
that of the star, q = 1 — M env /M e j. Clearly, the initial distribution of density in the progenitor will 
affect the final distribution of ejecta in the ejecta. 

There is one other parameter that plays an important role in determining the ejecta distri- 
butions from the convective outer envelopes of red supergiants: M env /M e j, the ratio of the outer 
envelope's mass to the total ejected mass. This ratio determines how much material (beneath 
the outer envelope) participates in the explosion. Although the outer envelope will wind up with 
the lion's share of the final kinetic energy, the mass of underlying ejected mantle is important in 
the blastwave stage. At the beginning of this phase, the explosion energy is contained within the 
mantle. During the course of the blastwave, as the explosion energy gets transferred to the outer 
envelope, the mass of the swept-up mantle material continues to affect the forward shock velocity 
(eq. 0). 

Note that the shock velocity in the outer envelope (as approximated by equation Q) is affected 
by the mass of the mantle, but not (to the accuracy of this approximation) by its structure. 
Indeed, there is reason to believe that this property should hold for the final ejecta distribution 
as well. Our pressure-based model for the ejecta is a local theory, making reference to the initial 
variables (m, ro,po) i n the progenitor. The distinction between remnant and mantle affects the 
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mass coordinate m, so the mass ratio M env /M e j affects our prediction for the outer envelope ejecta. 
But, the initial radius and structure of the mantle do not. So, we are justified in using simplified 
models whose mantles have a realistic mass but an idealized structure. 

The two parameters q = 1 — M env /M^ and M cnv /M e j are clearly correlated for real progenitors, 
because M± and M e j only differ by the remnant mass M rem . Indeed, only certain values of these 
two parameters are realistic, as depicted by Figure ||. Although the loss of the outer envelope mass 
in isolated stars and binaries is uncertain and intrinsically variable, the evolution of stellar cores 
and mantles appears to be nearly independent of mass loss and initial metallicity while a hydrogen 
envelope persists ( |Woosley fe Weaver 1995 , Maeder 1992| ). Assuming that the fiducial estimate 



of the remnant mass is correct, this implies that M rem and M man are well known as a function of 
initial stellar mass. The final value of the hydrogen envelope mass M env depends on uncertain mass 
loss, but it is unlikely that the final stellar mass M* of an RSG will exceed the initial main-sequence 
mass, except perhaps in some massive binaries. The circles in Figure |6] demonstrate the locus of 
the two parameters for stars that lose no mass; the stars on this line are the models of Maeder 
(1992), which also agree with those of Woosley &: Weaver (1995) and show no significant variation 
with metallicity. Varying M env between this maximum amount and zero, a star with a certain ratio 
M rem /M man will trace out a curve (shown as dotted lines) in the space of these parameters. As 
discussed by Woosley &: Weaver (1995), the mass loss rate of isolated stars increases with metallicity. 
Their calculated final states of solar-metallicity stars are shown as stars in the figure. 

It is evident from this figure that only a narrow swath of the parameter space of q and M env /M e j 
is likely to be inhabited by type II supernova progenitors. So, we are justified in focusing our 
attention on the one-parameter family (dashed line in Figure ^), 

^ = 1.12-1.08,, (48) 

for values of q between 0.3 and 0.8. For comparison, we have also considered values of M env /M e j 
greater by 0.05 at a given q, in order to capture the variation of quantities over the physical range. 

Our model stars have a dense central region that represents the mantle, which has a small 
radius (typically 10" 4 -??*), a uniform density, and the appropriate mass as given by M man /M e j. 
Above this lies a polytrope of index n = 3/2 that represents the stellar envelope; its structure is 
given by the value of q under consideration. Explosions in these progenitors are then simulated 
numerically as outlined in §||, and the results are analyzed to find the variation of the parameter a 
for the model of section |as a function of q. 

We base our trust in these models on demonstrations that their outer envelope ejecta are insen- 
sitive to the initial mantle distribution and can match the results of more sophisticated progenitor 
models. In Figure ||, we show that varying the mantle radius (and thus its density) alters the final 
distribution of mantle material slightly, but not the final distribution of the outer envelope material. 
This figure also shows that our models can match the outer envelope ejecta in the explosion of the 
15M model of Woosley, implying that large differences in the initial mantle distribution lead to 
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negligible differences in the final outer envelope ejecta. 

9. Calibration of models against simulations 

9.1. Simulation technique 

Our simulations are performed with an explicit, second-order, Lagrangian, finite-differencing 
code, as described by Richtmeyer h Morton (1967). Our code was written by Kelly Truelove, 
modeled after a code by John Holliman, as described in Truelove & McKee (1998). We have 
modified it slightly for the simulation of supernovae. Our progenitor distributions were kindly 
supplied by Ken'ichi Nomoto and Stan Woosley, and have been interpolated onto finer Lagrangian 
grids for more accurate dynamical results. For the study of lone polytropic envelopes, we solved 
the Lane-Emden equation numerically. After finding the solution corresponding to a desired value 
of q for a given value of the polytropic index n, we sampled the resulting structure at initial radii 
ro(i) oc i to generate the i'th Lagrangian element. The sampling parameter I was varied to capture 
the full range of mass, initial density and initial radius, and to give the fastest convergence for studies 
in which the number of zones was varied. When investigating the highest-velocity ejecta, we were 
motivated by spherical truncation of self-similar acceleration (see the Appendix) to distribute the 

1/3 

Lagrangian zones evenly in the parameter x Q . Runs were performed with between 500 and 10, 000 
zones, depending on the convergence properties and the fidelity desired for a particular result. 

Our calculations were performed with a maximum Courant number of 0.5, and with a coeffi- 
cient of the quadratic artificial viscosity law Cq = 2.5. These values were chosen to optimize the 
energy and entropy conservation, as well as the numerical shock structure. Small spikes in the 
ejected mantle of the BSG (seen in figure ||) are due to features (deviations from smoothness) in 
the low-resolution progenitor that we interpolated for our simulations; these are not caused by 
under-resolution of the hydrodynamics. Energy was conserved to within 3% in all our runs, with 
the discrepancy growing during shock propagation and holding steady after breakout. All calcula- 
tions were carried out until the thermal pressure was three orders of magnitude smaller than the 
kinetic energy density in every zone; at this stage, pressure gradients could safely be ignored. To 
determine the forward shock velocity we compared the entropy distribution with the initial den- 
sity distribution. This comparison was made before the passage of a reverse shock, for zones that 
experienced one. A slight loss of entropy conservation (decreasing with increasing resolution) was 
observed in the range of 10 — 100 breakout times, and we were careful to record the shock-deposited 
entropy of all zones by the time of breakout. 
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9.2. Comparisons against simulations of realistic progenitors 

We have compared in detail both the pressure-based and the harmonic-mean models for the 
final density distribution pft 3 against the output of the hydrodynamic code for the case of two 
progenitor distributions: a 15M red supergiant progenitor (Woosley & Weaver 1995, courtesy of 
S. E. Woosley), and a 16M Q progenitor for supernova 1987A (Shigeyama & Nomoto 1990, courtesy 
of K. Nomoto). The comparison is made graphically in Figures [8L [9| and 10, and the results are 
summarized in Table |(| 

Both models are typically more accurate in the hydrogen envelope of the red supergiant than 
that of the blue supergiant. This is to be expected, because the parameters of these models were 
tuned to match the results of simulation using polytropic analogs for convective stars (§0). However, 
the effect of the superadiabatic gradient on the high-velocity ejecta of the RSG is not included in 
these polytropic analogs. The pressure-based model, which includes more physics, better reflects 
the initial distribution of this material and only incurs ~ 34% error (weighted by log(t>/)) as a 
result. The harmonic-mean model, on the other hand, differs from simulation by about a factor 
of ~ 2.5 with the same measure. The difference can also be seen in the figures. This should be 
considered an error of the polytropic analog, or a demonstration of the effect of the superadiabatic 
gradient on the high- velocity ejecta. 

In the blue supergiant as well, variations of the effective polytropic index dominate the (log 
velocity-weighted) errors in the high-velocity ejecta. An inspection of the progenitor reveals that 
n e g approaches and then exceeds 3 in the outer half of the radius, and then becomes approximately 
constant in the outermost 10%. The simulation and the models both attain the desired self-similar 
power law index at the highest velocities. The deviations away from this index are due to both 
the effects of sphericity (§ |5.2| ) and the variation in polytropic index: steeper progenitors (higher 
^eff) gi ye shallower ejecta (lower l p 2) in self-similar theory (e.g., Chevalier &; Soker 1989| ). The 



pressure-based model does become shallower for lower n e fj because of the change in shock velocity, 
but not as shallow as the true numerical behavior. As a result, an error accumulates because the 
shape of the ejecta is slightly misrepresented, especially for material that originated in the outer 
quarter or so of the radius. 

It would be possible to expand the analytical description to include variations in n e g from 
their ideal values of 3/2 and 3 for red and blue supergiants. Such an analysis would be necessary 
for the association of observational features in the high- velocity ejecta with initial features of the 
progenitor. 

The mantle distribution of the red supergiant is badly represented in both models. The 
pressure-based model suffers from a gross underestimation of the entropy in this region, because 
we do not account for the effect of the reverse shock. The harmonic-mean model does not cover the 
RSG mantle at all, because it was fit only to the outer ejecta in analogous polytropic simulations. 



For the blue supergiant, the effect of the reverse shock is not as severe, and the pressure-based 
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model fares quite well. Because it is accurate in predicting the density of the mantle ejecta, the 
velocity of the hydrogen-helium boundary is also accurately predicted. The harmonic-mean model 
is aided by the weakness of the reverse shock as well, but this model includes no information about 
the difference between mantle and outer envelope. As a result, it incurs a significant error in the 
BSG mantle material. 

Note that the pressure-based model correctly predicts the density jump (eq. Jl5[]) between 
the hydrogen and helium for both red and blue supergiants. The effect of the reverse shock is 
apparent inward in mass from this jump; it depresses the central density of the red supergiant 
ejecta, producing an offset density maximum. We will discuss this phenomenon further in the next 
section. 

Besides the stellar models presented so far, we also investigated the explosion of model 3H11 
of Nomoto, Iwamoto & Suzuki 1995, provided to us by Ken'ichi Nomoto. This model, meant 
to describe the progenitor structure of SN193J, has an H envelope of mass O.IIMq, very small 
compared to the 3.4M ejected. The envelope's density structure is non-polytropic; indeed, there 
is a density minimum at the base of the envelope and a local maximum within it (seen in figure 6 of 
Nomoto, Iwamoto & Suzuki 195). The density then declines towards the surface, but shallows out 
again, with dlog p/dlog depth decreasing from 1 to 0.2. The reverse shock into the He mantle stalls 
after it has traversed a mass comparable to the envelope mass, leaving behind a jump in density. 
The shallow initial density law leads to steeply declining final density in the outermost hydrogen 
ejecta, as predicted by equation (p5|) and depicted by Suzuki & Nomoto (1995). The density jump 
and the variations of outer power law were not reproduced by our pressure-based or harmonic-mean 
models, although the helium mantle ejecta was fairly well represented. 



10. Conclusions 

In this paper we develop analytical approximations that capture the essential elements of 
spherically symmetric, adiabatic explosions in media of finite mass and extremely stratified density 
distributions: the progenitors of supernovae. Analytical predictions are motivated and then cal- 
ibrated by comparison to numerical simulations for realistic red and blue supergiant progenitors. 
In separate sections we address the velocity law of supernova shock waves (§|j), the expansion of 
the highest- velocity ejecta (§pl), and two types of models for the final distribution of ejecta. One 
of these models (§^j) specifies the pressure distribution in the ejecta, and derives the final density 
distribution; the other (§0) specifies the density directly as a function of final velocity. The benefits 
and drawbacks of these models, in application to red and blue supergiants, are summarized in Table 

0- 

We emphasize the structure of the outer (typically hydrogen) envelope ejecta in particular, 
because this layer is separated from the underlying mantle material by a large drop in initial 
density. This drop in density leads to an acceleration of the forward shock, followed by a reverse 
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shock into the mantle material as the forward shock decelerates. Although our models can account 
for the density jump in the ejecta caused by the acceleration of the forward shock (§ |6.1| ), we make 
no attempt to model reverse shock propagation. So, our models do not describe the mantle material 
accurately, although this drawback is much less severe for blue supergiants than for red ones. 

No doubt the most useful conclusion of this work is that the ejecta distributions from red 
supergiants are described by the single parameter q = 1 — M env /M+. This simplification occurs 
because these progenitors follow a sequence in parameter space (the dashed line in Figure |6|) as 
they lose mass, regardless of the mechanism of stellar mass loss. When these stars explode, their 
ejected hydrogen envelopes assume a distribution that can be fit reliably by the simple, harmonic- 
mean models of §0. The sequence of progenitors produces a sequence of ejecta distributions (with 
fits tabulated in Table ||). The simplicity of these harmonic-mean models makes them ideal for 
a systematic analysis of type II supernova light curves from RSGs. Such an analysis would be 
useful in calibrating the Expanding Photospheres Method, as discussed by Eastman, Schmidt and 
Kirshner (1996). 

However, this work is limited by the assumption of strict spherical symmetry in the analytical 
models and in the numerical simulations against which they are compared. The oblateness of a 
rotating progenitor (Chevalier & Soker 1989), the instability of accelerating shock waves (albeit 
slowly growing; Chevalier 1990, Luo & Chevalier 1994) and the Rayleigh- Taylor instability (e.g., 
Herant & Woosley 1994) have interesting effects that are not accounted for here. Moreover, by 
approximating supernovae as adiabatic, non-gravitating point explosions we have neglected physics 
that is important for the inner mantle material (where reverse shocks also limit our accuracy). 



Our model for the shock velocity (eq. fl7| ) is the simplest possible formula that combines 
powers of the local variables m, ro and po i n the progenitor and the blastwave energy E[ n in a 
manner consistent with both spherical and planar scaling laws, without introducing dimensional 
scales. The success (Figure |J) of this formula implies that blastwave shocks are controlled primarily 
by these scalings, even in situations where the postshock flow is quite complicated, including, at 
times, reverse shocks. This is very convenient, because it allows an accurate estimation of the 
entropy distribution deposited by this forward shock (although the entropy from a reverse shock 
is not accounted for). An accurate prediction of the shock velocity and its dependence on the 
stellar structure also permits a systematic analysis of the dynamics of the high-velocity ejecta 
(§H), including the maximum velocity achieved by the ejecta (§ l5.3| ) and an analysis of the burst of 
radiation that occurs at shock breakout (§ 5.3.3] ). 



We find in § 5.3.2 that a star must be sufficiently compact in order to produce any relativistic 
ejecta; for such stars, we make a rough estimate of the rest mass of relativistic ejecta. Applying 
these considerations to the supernova 1998bw and the gamma-ray burst, GRB 9804025, that may be 
associated with it (Soffita et al. 1998), we find that an explosion of relatively small ejected masses 
(IMq) and large energy (^ 10 52 ergs) could provide enough energy in relativistic ejecta both for 
the gamma-ray burst (Galama et al. 1998), and to power the inferred relativistically expanding 
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synchrotron shell (Kulkarni et al. 1998). Woosley et al. (1998) considered models for SN1998bw 
with larger ejected masses in order to fit the supernova light curve, and found accordingly smaller 
energies in the relativistic ejecta. Even assuming the existence of this ejecta, liberating its energy 
in a brief period requires an unusually high circumstellar density, as Woosley et al. point out. 

The velocity at which stellar material expands into space is related to the velocity of the shock 
front that set it into motion. However, the nature of this relationship is very different for material 
deep in the ejected envelope than for material near the stellar surface. Simple considerations (§|5|) 
reveal that the highest-velocity (initially outermost) material is cast away at some multiple of the 
shock velocity that traversed it, although the value of this factor varies with initial depth. At lower 
velocities, the smoothness of the pressure distribution controls the density distribution of the ejecta 
from deep within the progenitor. Here the shock velocity determines the final velocity indirectly, 
by means of the entropy distribution. The transition between these two behaviors is the result 
of a complicated, spherical rarefaction flow. But, simple models for the pressure distribution (§|6|) 
reproduce this transition quite accurately, especially in polytropic progenitors for which the high- 
velocity behavior is well understood (§ |5.2| ). Features in the ejecta at low velocities, most notably the 
structure of the ejected mantle, can be understood on the basis of approximate pressure equilibrium. 

One consequence of pressure equilibrium for the low-velocity ejecta is the presence a sharp 
density jump at the base of the outer envelope: the helium ejecta is denser and has a lower entropy 
than the hydrogen ejecta that surrounds it. This jump, quantified in § |6.l| , is greater in the ejecta 
from red supergiants (about a factor of 30) than from blue supergiants (about a factor of three). 
It results directly from the density difference between envelopes in a progenitor, and the shock 
acceleration that occurs there. As we have discussed, this feature is unlikely to be affected by 
the development of reverse shocks or the Rayleigh- Taylor instability, and may have interesting 
observational implications in the epoch of recombination. 

The much larger drop in density from helium to hydrogen between red and blue supergiants 
also leads to a much stronger reverse shock in RSGs than in BSGs. This limits the accuracy 
of our pressure-based models, which do not include the entropy deposited in the helium mantles 
(particularly in RSG mantles) by the reverse shock. It also has important consequences for the 
velocity of radioactive material in the ejecta. Consider, for instance, the fate of the innermost Q.bM & 
of ejecta in the 15M Q RSG and the 16M Q BSG progenitors depicted in Figures [7| and ||. In the 
numerical simulations, which include reverse shocks, this half solar mass is one to three times more 
dense at a given time in the RSG ejecta than in the BSG ejecta (the comparison is straightforward, 
since the ejected masses are quite similar). Because the outer envelope harbors almost all of the 
energy in the final state, its expansion velocity is essentially fixed. So, the expansion velocity of the 
mantle is determined by its density relative to the outer envelope, which depends on the difference 
in entropy between the two regions. This, in turn, is affected by the reverse shock. Our pressure- 
based models indicate what might happen if the reverse shock never developed, and the mantle 
were decelerated by a compression wave instead: it would then lack the reverse shock entropy, 
and be much denser (and expanding slower) in comparison to the outer envelope. Although these 
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concepts are the products of spherically-symmetric models, similar considerations should apply in 
realistic explosions as well. 

It is noteworthy that the reverse shock develops inward (in mass) from the density jump 
between the helium and hydrogen ejecta, so that there is a density peak at the outer helium ejecta 
in red supergiants (Figure |?j) . 

In BSGs the reverse shock is much less severe than for RSGs, and (at least for the model 
16Mq BSG provided by Nomoto) seems to develop far inward in mass from the hydrogen-helium 
boundary. The pressure-based model tracks the mantle ejecta quite well. Indeed, the difference 
between the mantle and the outer envelope can be ignored entirely without too much error: such is 
the spirit of the Chevalier & Soker (1989) model for SN1987A, whose elements we have incorporated 
into our harmonic-mean model for BSGs. Here, the mantle and the inner region of the hydrogen 
ejecta can be grouped together as a single region of gently declining density. 

Regardless of the mantle, it appears that the ejecta density distribution in the bulk of the 
hydrogen envelope is a flat or gently declining function of velocity (or mass) for both red and blue 
supergiants. At some break velocity, the structure of Pf(vf) rolls over toward the steeply declining 
power law it must assume at high velocities. For red supergiants, a quantitative statement can be 
made: density falls by only a factor of 20 — 40% in the inner half of the hydrogen ejecta mass. (This 
density drop decreases with increasing mass ratio q, despite the fact that the inner power law l\ 
gets steeper [Table but note that the hydrogen envelope is a smaller fraction of the ejecta as q 
increases.) 

Chevalier (1976), making an analogy with Sedov blastwaves of different indices, suggested that 
a detached shell (density maximum) should form in the ejecta at shells for which din po/dlnro 
decreases smoothly through the value —2.53 in the progenitor. If this were strictly true, a detached 
shell would form in the outer envelope ejecta of all the progenitors we have considered. No such 
density maxima actually formed in the outer envelopes, although density maxima did occur in the 
ejected mantles (of RSGs especially, as described above). The presence of density maxima in the 
ejected mantles depends on the strength of the reverse shock, and does not appear to be related to 
Chevalier's proposed mechanism. 

It is clear from our simulations that when a progenitor deviates from a polytropic structure 
in the outer half of its radius, perhaps because of inefficient convection or the details of radiative 
transfer, the resulting high-velocity ejecta distribution bears the imprint of this effect. This is 
apparent when the ejecta from simulations with realistic stellar progenitors are compared against 
the ejecta of polytropic models. Our analysis of the high-velocity ejecta of polytropic models is 
limited in applicability by this effect. On the other hand, the relative constancy of the shock 
velocity parameter (3\, the postshock acceleration ratio (vf/v s ) p , and the truncation radius Rt for 
planar acceleration all suggest that it may be possible to extend our semi-analytical theory to 
accommodate variations in n e fj. With such a theory, it may even be possible to constrain models 
of inefficient convection with the early light curves of type II supernovae. 
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Our pressure-based model, being derived from a simple model for the entropy deposited in the 
blastwave shock, is a local theory that relates the final density pj of a mass shell with its initial 
density po, initial radius rg, and enclosed ejected mass m. A strong reverse shock, whose passage 
cannot be predicted by these local variables alone, invalidates this model in the mantles of stars. 
But, the local theory is exceptionally accurate where it is valid. 

Our harmonic-mean models, on the other hand, are global models for the ejecta distributions 
of supernovae, in the sense that they specify the general run of final density versus velocity that 
results from a specific type of progenitor structure. These models are appealing for their simplicity, 
although they contain little information about the dynamics of envelope expulsion in supernovae. 

The authors are indebted to Stan Woosley and Ken'ichi Nomoto for providing supernova 
progenitor structures, and to J. Kelly Truelove for providing the computer code that we modified 
for the purposes of this paper. CDM thanks John Monnier, Ben Metcalf and George McLean for 
many valuable suggestions, Lars Bildsten for clarifying stellar structure and evolution and reviewing 
a draft of this work, and Douglas Leonard for discussing light curves. We also wish to thank our 
referee, Roger Chevalier, for insightful and helpful comments. CDM was supported in part by a 
National Science Foundation Graduate Fellowship and an ARCS fellowship. The research of both 
CDM and CFM is supported in part by the National Science Foundation through NSF grant AST 
95-30480. 

A. Fully Lagrangian, self-similar treatment of shock emergence 

We present a calculation of the self-similar, planar flow in a power law density distribution with 
an edge, the problem posed by Gandel'man & Frank-Kamenetsky (1956) and solved by Sakurai 
(1960). Unlike previous treatments, we use a Lagrangian approach for the pre-breakout flow as 
well as the post-breakout flow. This has the advantage that variables can be compared against 
their postshock values, rather than against their breakout values (which previously were related 
to postshock values with an Eulerian treatment). To do this, we first define the space and time 
coordinates x = 1 — r/R+ and At = t — t (breakout). The time at which a mass shell experiences 
the shock is Ato(m). For our self-similar coordinate we choose rj = At/ Ato(m), and for the flow 
variable we choose S(tj) = x(m)/xo(m). Our variables rj and S(rj) are both unity in the initial 
state, and decrease toward — oo as the flow progresses. Breakout occurs when At = rj = 0. 

To maintain consistency with previous authors, we define the initial density distribution as 
Po(m) = fci[i?*xo(m)] n , and the shock velocity to be v s = —k2{R*xo(m)]~ x , where A = n(3\. Then, 
At = -(R+x ) 1 - x /k 2 (l + \). 

Considering the time variable At to be parameterized by r\ and xo, the constraint of constant 
time, dAt = 0, gives drj = —r](l — X)dxo/xQ. This relation allows us to take the necessary derivatives 
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to solve for the flow. First of these is the density, defined in planar geometry by: 

P (V, *o) = Po (§) m = m - J ^ m . (Al) 

Now, we evaluate the pressure as p(r], xq) = p s (p/p s )' y , where p s and p s are the postshock pressure 
and density. The acceleration in planar geometry is, R+x = — {dp/ dr)j\ t / P = — (dp/drQ) / \ t /po- From 
the definitions of S(t]) and rj, on the other hand, x = xqS" (77) /t§ . Combining the two expressions for 
the acceleration, we have the following differential equation for S(r]) to describe planar, self-similar 
flow: 

, 2c „ 2 fj-iy ( n-2X ^ \S' + {X-l)r]S" 



W " = — [—1) I p-tA+W ^-^ p-tA + l^ l- (A2) 

The initial condition for this equation is at the time 77 = 1, when a mass shell is hit by the 
shock wave. Because the mass shell is in its initial position at this time, and because its postshock 
density and velocity are given by the Rankine-Hugoniot jump conditions for strong shocks, 

"M- 1 - s ' (1) " hTWTT)- (A3) 

Note that the ratio v/v s at any point in the flow is equal to (A + l)S'(rj), for the value of 77 
appropriate to the time and mass shell considered. 

We used the technique described by Sakurai to determine A (or equivalently j3\) for 7 = 4/3 
and n = 3/2 and 3. Then we solved for the evolution of a mass shell's position using equation ( |A2[ ), 
starting with the initial conditions of equation ( |A3| ). We continued the numerical solution out to 
\rj\ > 10 8 , and used the known asymptotic solution to find the limit S'(r] — > —00). For the purpose 
of studying the solution, we then fit the resulting velocity history v/v s in powers of (2 — r?) -1 / 3 . 
This parameterization was chosen to show the asymptotic dependence of S{rj) as rj — > —00, and 
also to be finite for all 77 < 1. We obtain, 

v(m,ri) 1.0656 0.5872 1.1235 . 

V ~ 2.1649 - -tttt + \~<77q (n = 3/2 , 

v s (m) (2-r ? ) 1 /3 (2-r/) 2 /3 2 - r] y ' h 

v(m,ri) 0.7638 0.2535 0.7861 . . , AJ . 

v / '' ~ 2.0351 - -ttt + -77^ (n = 3). (A4) 

v s (m) (2-r?) 1 /3 (2-r?) 2 /3 2 - r? v ; v ' 

These fits, meant to capture the late-time evolution accurately, are accurate to within 2% and 
1.5%, respectively, for r\ < —1, although they miss the breakout value (vf/v s = 6/7) by 34% and 
14%, respectively. For the purpose of studying the truncation of the self-similar acceleration law, 
as discussed in §|| it is useful to make similar fits based on the normalized distance coordinate 
S(rj) = x/xq rather than the time coordinate 77: 

v(m,S) A 1.2504 0.8473 0.9817 , 

< m ^ S ) „ 2 0351 °- 826Q I °- 3 ° 94 °- 5244 fn-31 fA5) 
v s (m) 2 '° 351 (2 -5)V3 + ( 2 -5)2/3 2-S {n - d} - (A5) 



-39- 



Much like the previous fits, these are both valid to within 1% for S < — 1, but they miss the 
breakout value by 9% and 16%, respectively. 

Litvinova & Nadezhin (1990) suggested that spherical effects would act to truncate the planar 
acceleration once a shell had reached a truncation radius Rt ^ i?*. Using our solution to the 
self-similar expansion, we can evaluate this prescription for different truncation radii and compare 
it to the results of fully spherical calculations (§[5^). The normalized position St corresponding to 
R t is, 

S t = - fr/**- 1 (A6) 

X 

Taking Rt = 3R± as an illustrative example (Figure |3|), the self-similar solution can be fit as a 

1 /Q 

polynomial in x in accordance with equation ( |A6| ): 

- (—) fl- 0.4522xn /3 + 0.1717xn /3 - 0.04328x0) (n = 3/2), 

v s {ro) \v s Jp\ u u / 

^T^l - (—) fl - 0.3174xn /3 + 0.0373xn /3 + 0.0199x0) (n = 3). (A7) 

v s {r ) \v s J p \ > 

These fit the self-similar solution quite well, with maximum errors 0.3% and 0.5% for n = 3/2 and 
3, respectively, in the entire range < xo < 1. Clearly these are quite similar to the outcome of 
full spherical dynamics as given by equation (^). 
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mass (solar units) 

Fig. 1. — The red and blue supergiant progenitors used in this paper. The 15M RSG (solid line, 
Woosley & Weaver 1995) was provided by Stan Woosley. The 16M & BSG (dashed line, Shigeyama 
&; Nomoto 1990) was provided by Ken'ichi Nomoto. The boundaries between remnant and ejecta 
(dash-dot line), and between mantle and outer envelope material (M a , circles), are marked in both 
cases. We take the remnant to be the edge of the iron core; this gives remnant masses of 1.48 and 
1.69M Q for the RSG and BSG, respectively. 
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(m-M ) /M 

rem e j 

Fig. 2. — The forward shock velocity in a numerical simulation of the 15Mq RSG progenitor S15A 
of Woosley & Weaver 1995 (solid line), and a comparison between our shock velocity formula (eq. 
]i~7| , dotted line) and the formulae of Kompaneets (1960, thick dashed line), Klimishin & Gnatyk 
(1982, thin dashed line), and Koo & McKee (1990, dash-dot line). Except for our formula, shock 
velocities have been normalized to match the numerical value at rh = 0.7. The constant R s j in the 
model of Koo & McKee has been chosen to optimize the shock velocity in the outer envelope. To 
be fair, ours is the only formula designed specifically for the context of super novae. 
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Fig. 3. — The ratio of shock to final velocity for a sequence of simulated explosions in polytropes 

1/3 

without mantles (index n = 3/2, q = 0.02, 0.3 and 1), parameterized by the radial variable x . If 
the progenitors were homologous copies of one another, dimensionless quantities such as this would 
not vary. The progenitors are most similar at the outermost radii, a feature that holds for this 
ratio as well. The ratio in each simulation approaches the planar, self-similar value (dash-dot line). 
The suggestion of Litvinova & Nadezhin (1990) that spherical expansion may truncate self-similar 
expansion at a radius E4 gives a definite prediction of the variation of this ratio with depth (dashed 
lines, for R t = 1.5, 2 and 3R*). Our fit to the dependence of Vf/v s with x$ is also plotted (eq. fl26| , 
dotted line). 
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Fig. 4. — Similarity of the outer regions of polytropic envelopes is reproduced as similarity of the 
high- velocity ejecta. a: The sequence of polytropic progenitors used to simulate red supergiants 
(§||), plotted as normalized initial density po/pi against the variable ro/(R± — ro), in which po/pi is 
a power law for the region in which polytropic envelopes share a common form (§|3.1|). b: The ejecta 
distribution produced in these models: normalized final density against normalized final velocity. 
For comparison, the self-similar planar power law is also shown (dashed line). 
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Fig. 5. — This figure demonstrates the applicability of polytropic models to red supergiant en- 
velopes, as well as the insensitivity of the outer envelope ejecta to the size and structure of the 
mantle. The red supergiant model of Woosley & Weaver (1995, solid lines) has a mantle radius 
i? man ~ 0.02i?^ and a complicated mantle structure resembling po oc r~ 3 . The polytropic models 
have mantles of constant density, with R man = 10 _3 i2* (dashed lines) and i? man = 10~ 4 i?* (dotted 
lines). Nevertheless, the outer envelope ejecta is well-modeled by both progenitors. The superadia- 
batic gradient is significant in the subsurface layer of the progenitor, which produces a high-velocity 
density distribution that differs from polytropic models. 



-48- 




1-M /M = q 

env * 

Fig. 6. — The mass ratios that parameterize type II supernovae. Plotted is the envelope fraction 
in the ejecta, versus the ratio q of core to stellar mass. If remnant masses are correctly estimated, 
progenitors must lie in the swath between the solid lines regardless of mass loss. As a star loses its 
envelope, it travels from the point corresponding to no mass loss (circles, from Maeder 1992, for 9, 
12, 20, 25, 40 and 60 M & ) along a curve (solid lines, for 9 and 4OM ) towards the point (1, 0). The 
results of lone stellar mass loss at solar metallicity (stars, Woosley & Weaver 1995) for 9, 12, 20, 
and 25 Mq lie in the interior. Stars initially above about 27Mq will only retain their envelopes if 
the metallicity is low. The remnant mass cannot be negative (hatched region), but it could be zero 
in a thermonuclear explosion. Dashed line: the sequence of models (eq. ^8|) considered in §||. 
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0.5 1 0.5 1 5 10 

m/M v /v 

e j f * 

Fig. 7. — A comparison of the pressure-based model (§|6|, dotted lines) with the result of simulation 
for realistic progenitors (solid lines, with circles at the H-He boundary). The maximum ejecta 
velocity is marked (§ |5.3| , dashed lines). Left side: normalized final density pft 3 /p*t^, with respect 
to the normalized ejecta mass rh = m/M e y Right side: ptt s /p±t^ versus the normalized velocity 
Vf/v±. The 15Mq red supergiant progenitor model of Woosley & Weaver (1995) was the input 
into both the simulation and the analytical model. The parameter a was chosen according to the 
formulae in Tables |I| and ||. Note that the models fail interior to the base of the hydrogen envelope, 
because of the entropy deposited in the reverse shock, but the jump in density at the base of the 
outer envelope is correctly estimated (§|6.1|). 
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m/M v /v 

e j f 



Fig. 8. — A comparison of the pressure-based model with simulation, but here the progenitor was 
taken to be the 16M© blue supergiant model of Shigeyama & Nomoto (1990). As in Figure 0, the 
solid line is the simulation, with circles at the He-H boundary; the dotted line is the output of 
the pressure-based model, and a dashed line marks the maximum velocity to which the adiabatic 
approximation is valid, which is approximately the speed limit for the ejecta. Note the lesser effect 
of the reverse shock for the BSG compared to the RSG, the higher maximum velocity, and the 
greater range of ejecta densities. Again, the jump in density at the H-He boundary is correctly 
determined by the model. 
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e j f 



Fig. 9. — Application of the harmonic-mean model (§|7] and Table [|: dotted lines) to numerical 
simulation (solid lines, with circles at the H-He boundary), for the RSG progenitor of Woosley 
&i Weaver (1995). Again, the dashed line shows the maximum velocity of the ejecta. For RSG 
progenitors, the model is derived from a polytropic progenitor, as tabulated in Table It is 
meant to describe the hydrogen envelope only. Since the models prescribe Pf(vf), the model mass 
coordinate is obtained by integrating inward from rh = 1 at high velocities. 
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Fig. 10. — A comparison of the harmonic-mean model (§^: dotted lines) with simulation (solid 
lines), here for the blue supergiant progenitor of Shigeyama & Nomoto (1990). For the BSG 
progenitor, the constraints of mass and energy conservation are imposed, and the inner power-law 
is chosen to agree with the self-similar rarefaction model of Chevalier & Soker 1989 (Table |5|). This 
model is therefore meant to cover the entire ejecta distribution. 
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Table 1. Models for the ejecta density. 





Pressure-Based Models (§||) 


Harmonic-Mean Models (§^) 


Model output pf(m)t 3 


p f (v)t 3 


Type of model Local, Lagrangian 


Global, Eulerian 


Flexibility Responds to detailed stellar structure Fixed analytical form 


Benefits 


Gives jump in pj at mantle 


Simple 


Drawbacks Requires initial structure 


Insensitive to detailed structure 


RSG mo 


dels^ Parameter a fit as a(q) 


3 parameters fit with q 




H envelope most accurate 


H envelope only 


BSG, lb, 


Ic models"'" Approximate a = 


Generalized model of CS89* 




Applies to mantle 


Approximates mantle 


^Models for red supergiants are based upon polytropic models for the progenitors; these are 


primarily determined by one parameter, q. 




^Radiative progenitors are not easily modeled with simple progenitors, but this is counterbalanced 


by the fact that features in their ejecta are less pronounced. 




*CS89: 


Chevalier & Soker 1989. 






Table 2. Quantities associated with self-similar, planar expansion. 


Quantity 


n = 3/2 n = 3 


Comments 


A 


0.2863 0.5573 


Method of Sakurai (1960) 


01 


0.1909 0.1858 


Method of Sakurai (1960) 




0.2071 0.2071 


Appx. of Whitham (1960) 


2.1649 2.0351 


Derived in Appendix 




3/5 3/4 


Initial density distribution 


Wv s 


-0.1145 -0.1394 


Equation (||) 


w s 


-0.3218 -0.3965 


Equation (||) 


Wp 


1.3436 1.4181 


Equation (||) 




1.3624 1.3624 


Equation (|23|) 
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Table 3. Formulae for high-velocity distributions of ejecta quantities. 



Quantity High-velocity formula Validity 



Convective outer envelopes (n = 3/2): 

Vf/v* 2.06(pi/^)-°- 08 (l - m)~ ail45 (l + 0.305m) 5m £ 0.8 

p f t 3 /p±tl 0.079(/>i/p*)-°- 23 (l-m) L344 (l + 0.675m) 6m £ 0.8 

P/i 4 /W 4 0.0027(pi/p*) - 02 (l - m) L362 (l + 0.965m) 5m g 0.7 

Radiative outer envelopes (n = 3): 

Vf/v* 1.92(p 1 /^)-° 05 (l - m)-°- 1394 (l - 0.235m - 0.0955m 2 ) 5m £ 0.9 

p f t 3 /p±tl 0.081(pi/p^)°- 14 (l-m) 1 - 418 (l + 0.495m + 0.615m 2 ) 5m ^ 0.9 

p f t 4 /p±ti 0.0026(/9i/^) - 01 (l-m) L362 (l + 1.035m + 0.375m 2 ) 5m £ 0.75 



Note. — Final distributions depend on the parameter pi/p* defined by equation (||), describing 
the subsurface layers of the progenitor (§ |3.l|) . Characteristic values such as v+ are combinations 
of the ejecta mass, energy, and radius (§§)• For each quantity, the distribution limits to a power 
law at the highest velocities. The factor that depends on 5m describes the deviation from this 
power law toward lower velocities, primarily due to the onset of spherical expansion. Here, 5m = 
[(1 - m)Goi//o*)- 1 ] 1 / 3 < n+1 >. These formulae are valid in the region for which envelopes of the 
same value of n, but different internal structures, share a common form in numerical simulations. 
Significant digits are apportioned according to accuracy: most accurately determined are the power 
law exponents, i.e., the exponents on (1 — m). 
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Table 4. Parameters for 


modeling the H ejecta 


from RSGs. 


P a V £1 TY1 P"i" PT* 


Value 


lypical KbG sequence (eq. Pq ): 






Q 


5.58g 2 -0.10g- 


1.17 


„ j.3 / ~ -z-3 


10~ 3 (5.21g 2 - 10.12(7 + 5.40) 


break/ ^* 


0.93<? + 1.94 






-3.06g 2 + 2.03g 


- 0.46 




-11.731 






4.5 




Pbreak^M** 


10" 4 (-9.87g 2 + 


4.67g + 4.87) 


break/ ^* 


0.80<? + 1.85 




^1 


1.68g - 1.89 






-11.895 






3.8 




Partial derivatives, for extrapolation: 






da/d(M env /M ei ) 


-17.9g + 3.67 




d{ PhrcSLk t 3 /p^)/d(M env /M ei ) 


10- 3 (-10.12g + 


9.53) 


9{v p brcak/^ )/d{M env /M e j ) 


1.14g - 1.43 




d/ p i/d(M cnv /M ej ) 


3.82g - 1.38 




a(pbreakt 4 M^)/9(M cnv /M ej ) 


10- 3 (1.53<? - 0.43) 


<9( V p break A>* )/d(M cnv / M ej ) 


0.28<? - 0.81 




dl pl /d(M env /M ei ) 


-0.74g-2.13 





Note. — The parameter a of equation (|4T|) , and the parameters used in equations (f46| ) and (^7|). 
The latter specify the harmonic-mean approximation for the hydrogen ejecta of red supergiants 
described in §[|; in the former, a specifies the more complicated pressure-based model of §^. Pa- 
rameters that depend on q are derived from the one-parameter family of progenitors described by 
equation ( p8| ) and shown as the dashed line in Figure |(| Partial derivatives corresponding to chang- 
ing the remnant mass without changing M+ or M env (i.e., at fixed q) are derived from comparison 
with a second sequence of models offset by 0.05 in M env /M e j. The fit for a is valid in the range 
0.3 < q < 0.6; others are valid for 0.3 < q < 0.8. 
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Table 5. Fiducial ejecta models for radiative progenitors. 



Parameter 


Value 


Comments 


a 





Typical value 


Pbrcak* 3 / P*tl 


1.91 x 1(T 3 


Energy and mass conservation 


Vp break/ V± 


2.30 


Energy and mass conservation 


Ipl 


-1.06 


Theory of Chevalier & Soker 1989 


lp2 


-10.176 


Self-similar value (n = 3) 


Vp 


~ 4.5 


From RSG models 


Pbrcaki 4 /W^ 


~ 2.4 x 10~ 4 


Gleaned from simulation 


Vp break/ V* 


~ 2.3 


From RSG models 


Ipl 


2.22 


Theory of Chevalier & Soker 1989 


lp2 


-9.77 


Self-similar value (n = 3) 


Vp 


3.8 


From RSG models 



Note. - - Parameters for modeling the entire ejecta of radiative stars. Without a sequence of 
reliable models to fit, these models are only educated guesses. The harmonic-mean model (all 
but a) is a combination of the self-similar rarefaction hypothesis of Chevalier & Soker 1989 with 
self-similar planar expansion for n = 3, and softening parameters y p , y p taken from the RSG fits. 
We simply assume a = to specify a pressure-based model. The pressure-based model reflects the 
detailed progenitor structure through the entropy deposited in the shock, whereas the harmonic- 
mean model is inflexible. Note that the effect of the reverse shock is less severe for these progenitors 
than for RSGs, so the mantle is better represented. 
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Table 6. Calibration of ejecta models for realistic progenitors. 



Error estimate 


RSG* 




BSGt 




H envelope 


Mantle 


H envelope 


Mantle 


Pressure-based model (§§) 










( 5 Pf/pf)rh 


0.07 


6.1 


0.10 


0.71 


[5pf(fh)/pf(m)] max 


2.1 


21 


7.5 


5.8 


( S Pf/pf)log(v f ) 


0.34 


4.5 


0.96 


3.0 


[8pf{vf)/pf{v f )\ max 


1.4 


32 


4.5 


15 


Harmonic-mean model (§[?]) 










( S Pf/Pf)rh 


0.20 




0.31 


0.70 


[6pf(fh)/pf(rh)] max 


4.5 




0.83 


22 


( 5 Pf/pf)lo g (v f ) 


1.5 


5.0 


0.53 


3.2 


[Spf(vf)/p f (vf)] max 


10 


18 


2.4 


19 



*15M RSG of Woosley & Weaver 1995, courtesy of Stan Woosley 
^16Mq BSG of Shigeyama & Nomoto 1990, courtesy of Ken'ichi Nomoto. 



Note. — Mean and maximum relative errors in the density models of this paper, for both the 
Lagrangian and Eulerian functions pf(m)t 3 and pf(vf)t 3 . Mean errors are computed according to 

the formula (Spf/pf) x = exp (^[log(p mo d e i(x)/ Pf(x))] 2 ^ — 1> where ( ) x is a mean weighted evenly 
in the variable x. The maximum error is defined consistently: an error of 1 implies a factor of 2 
(in either direction) between pf(x)t 3 and p m odei(x)t s . In RSGs, the harmonic mean approximation 
does not cover the mantle. Errors in the H envelope are dominated by variations in the effective 
polytropic index from its ideal value (3/2 for RSGs and 3 for BSGs) in the outer layers of the 
progenitor. Errors in the mantle are dominated by the effect of the reverse shock, which is not 
addressed in these models. 



